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1 .  ] n I  reduction 

Consider  the  evaluation  of  'ho  product  of  n-1  matrices 

M  =  M,  X  M,  X  •  X  M  ,  (1 ) 

12  n-1 

where  M.  is  a  w.  x  w.  ,  matrix.  Since  matrix  multiplication  satisfies 

iii  +  l 

the  associative  law,  the  final  result  M  in  (1)  is  the  same  for  all  orders 
of  multiplying  the  matrices.  However,  the  order  of  multiplication 
greatly  affects  the  total  number  of  operations  to  evaluate  M  .  The 
problem  is  to  find  an  optimum  order  of  multiplying  the  matrices  such 
that  the  total  number  of  operations  is  minimised.  Here,  we  assume  that 
the  number  of  operations  to  multiply  a  px  q  matrix  by  a  (]  Xr  matrix 
is  pq r  . 

In  [  1 3 [7 3  ,  a  dynamic  programming  algorithm  is  used  to  find  an 

3  2 

optimum  order.  The  algorithm  needs  0(n  )  time  and  0(n  )  space.  ln[2J, 
Chandra  proposed  a  heuristic  algorithm  to  find  an  order  of  computation  which 
requires  no  more  than  2T^  operations  where  T^  is  the  total  number  of  opera¬ 
tions  to  evaluate  (1)  in  an  optimum  order.  This  heuristic  algorithm  needs 
only  O(n)  time.  Chin  |^3j  proposed  an  improved  heuristic  algorithm  to  give  an 
order  of  computation  which  requires  no  more  than  1.  25  T^  .  This  improved 

heuristic  algorithm  also  needs  only  O(n)  time. 

* 

In  this  paper  we  first  transform  the  matrix  chain  product  problem  into 
a  problem  in  graph  theory  -  the  problem  of  partitioning  a  convex  polygon  into 
non-intersecting  triangles,  see  [9][  10][  1  l]t  12]  ,  then  we  state  several  theo¬ 
rems  about  the  optimum  partitioning  problem.  Based  on  these  theorems,  an 
0(n  log  n)  algorithm  for  finding  an  optimum  partition  is  developed. 


2 .  1  'a rtitioning  a  convex  polygon 

Given  an  n- sided  convex  polygon,  such  as  the  hexagon  shown  in 
Fig.  1,  the  number  of  ways  to  partition  the  polygon  into  (n- 2)  triangles 
by  non-intersecting  diagonals  is  the  Catalan  numbers  (see  for  example, 
Gould  [8]).  Thus,  there  are  2  ways  to  partition  a  convex  quadrilateral, 

5  ways  to  partition  a  convex  pentagon,  and  14  ways  to  partition  a  convex 
hexagon. 

Let  every  vertex  V.  of  the  polygon  have  a  positive  weight  w^  .  We 
can  define  the  cost  of  a  given  partition  as  follows;  The  cost  of  a  triangle 
is  the  product  of  the  weights  of  the  three  vertices,  and  the  cost  of  parti¬ 
tioning  a  polygon  is  the  sum  of  the  costs  of  all  its  triangles.  For  example 
the  cost  of  the  partition  of  the  hexagon  in  Fig.  ]  is 

w,w_w_  +  w.w-w.  +  w,w,w.  f  WWW,  . 

123  136  346  456 


(2) 


If  wc  erase  the  diagonal  from  to  and  replace  it  by  the  diagonal 

from  Vj  to  ,  then  the  cost  of  the  new  partition  will  he 

w.  w  w  +  w ,  w  w  t-  w,  ww,  4-  ww  w.  .  (3) 

123  134  146  456  '  ' 

We  will  prove  that  an  order  of'multi plying  (n-1)  matrices  corre¬ 
sponds  to  a  partition  of  a  convex  polygon  with  n  sides.  The  cost  of  the 
partition  is  the  total  number  of  operations  needed  in  multiplying  the 
matrices.  For  brevity,  we  shall  use  n-gon  to  mean  a  convex  polygon 
with  n  sides,  and  the  partition  of  an  n-gon  to  mean  the  partitioning  of  an 
n-gon  into  (n-2)  non- intersecting  triangles. 

For  any  n-gon,  one  side  of  the  n-gon  will  be  considered  to  be 
its  base,  and  will  usually  be  drawn  horizontally  at  the  bottom  such  as 
the  side  V  -V  in  Fig.  1.  This  side  will  be  called  the  base,  all  other 

i  6 

sides  are  considered  in  a  clockwise  way.  Thus,  V^-V^  ’s  t*ie  first 
side,  Vz-  V3  the  second  side,  .  . .  ,  and  V^-  the  fifth  side. 

The  first  side  represents  the  first  matrix  in  the  matrix  chain  and 

the  base  represents  the  final  result  M  in  (1).  The  dimensions  of  a  matrix 

are  the  two  weights  associated  with  the  two  end  vertices  of  (he  side.  Since 

the  adjacent  matrices  are  compatible,  the  dimensions  w^xw^,  w^Xw^, 

...,w  ,  X  w  can  be  written  inside  the  vertices  as  w, ,  w  w  .  The 

n-1  n  1  2  n 

diagonals  arc  the  partial  products.  A  partition  of  an  n-gon  corresponds 
to  an  alphabetic  tree  of  n-1  leaves  or  the  parenthesis  problem  of  n-1 
symbols  (see,  for  example,  Gardner  [6]).  It  is  easy  to  see  the  one-to- 
one  correspondence  between  the  multiplication  of  n-1  matrices  to  cither 


the  alphabetic  binary  tree  or  the  parenthesis  problem  of  n-1  symbols. 
Here,  we  establish  the  correspondence  between  the  matrix-chain  product 


and  the  partition  of  a  convex  polygon  directly. 

Lemma  1.  Any  order  of  multiplying  n-1  matrices  corresponds  to  a 
partition  of  an  n-gon. 

Proof.  We  shall  use  induction  on  the  number  of  matrices.  For  two 
matrices  of  dimensions  w^Xw^>  w2*  w3*  ^cre  *s  or|ly  one  way  of  multi¬ 
plication,  this  corresponds  to  a  triangle  where  no  further  partition  is 
required.  The  total  number  of  operations  in  multiplication  is  wjw2w3* 
the  product  of  the  three  weights  of  the  vertices.  The  resulting  matrix  has 
dimension  w^  X  w^  .  For  three  matrices,  the  two  orders  of  multiplication 
(MjXM^)X  and  correspond  to  the  two  ways  of  parti¬ 

tioning  a  4-gon.  Assume  that  this  lemma  is  true  for  k  matrices  where 


k  £  n-2  ,  and  we  now  consider  n-1  matrices.  The  n-gon  is  shown  in 


Let  the  order  of  multiplication  be  represented  by 

M  =  (M.  xM  X*  •  •  xM  . )  X  (M  x*  •  ■  XM  ) 
12  p- 1  p  n- 1 


i.e.,  the  final  matrix  is  obtained  by  multiplying  a  matrix  of  dimension 

(w.Xw  )  and  a  matrix  of  dimension  (w  Xw  ).  Then  in  the  partition  of  the 
1  p  P  n 

n-gon,  we  let  the  triangle  with  vertices  Vj  and  have  the  third  vertex  V  . 

The  polygon  V^-V^-  •  •  •  -V^  is  a  convex  polygon  of  p  sides  with  base  V^-V 

and  its  partition  corresponds  to  an  order  of  multiplying  matrices  M^,  . . . ,  M^  ^ , 

giving  a  matrix  of  dimension  w^Xw^.  Similarly,  the  partition  of  the  polygon 

V  -V  V  with  base  V  -V  corresponds  to  an  order  of  multiplying 

p  p+1  n  p  n 

matrices  M . M  ,,  giving  a  matrix  of  dimension  w  xw  .  Hence  the 

p  n-1  p  n 

triangle  V  V  V  with  base  V  -  V  represents  the  multiplication  of  the  two 
1  p  n  In 

partial  products,  giving  the  final  matrix  of  dimension 


w.  xw 


Lemma  2.  The  minimum  number  of  operations  to  evaluate  the  following 
matrix  chain  products  are  identical. 


M  ,  X  M  ,  X  •  •  •  X  M  ,  X  M 
1  2  n- 2  n-1 


M  X  M,  X 

n  1 


XM  ,  X  M 
n-  3  n 


2 


M  X  M ,  X  •  •  •  X  M  ,  X  M 
2  .3  n-1  n 

where  M.  lias  dimension  w.  X  w.  .  and  w  ,  ,  =  w,  .  Nolo  (hat  in  the 
i  l  t  H  n+ 1  1 

first  matrix  chain,  the  resulting  matrix  is  of  dimension  w^  hy  wn  •  'n 
the  last  matrix  chain,  the  resulting  matrix  is  of  dimension  w^  by  w^  . 

But  in  all  the  cases,  the  total  number  of  operations  in  the  optimum  orders 
of  multiplication  is  the  same. 


6 


Proof.  The  cyclic  permutations  of  the  n- 1  matrices  all  correspond  to 
the  same  n-gon  and  thus  have  Ihe  same  optimum  partitions.  ■ 

(This  Lemma  was  obtained  independently  in  f  4]  with  a  long  proof.  ) 

From  now  on,  we  shall  concentrate  only  on  the  partitioning 
problem. 

The  diagonals  inside  the  polygon  are  called  arcs.  Thus,  one 
easily  verifies  inductively  that  every  partition  consists  of  n-2  triangles 
formed  by  n-3  arcs  and  n  sides. 

In  a  partition  of  an  11-gon,  the  degree  of  a  vertex  is  the  number 
of  arcs  incident  on  the  vertex  plus  two  (since  there  are  two  sides 
incident  on  every  vertex). 


Lemma  3.  In  any  partition  of  an  n-gon,  n  2  4,  there  are  at  least  two 
triangles,  each  having  a  vertex  of  degree  two.  (For  example,  in  Fig.  1,  the 
triangle  ^as  vertex  ^2  degree  2  and  the  triangle  has 

vertex  with  degree  2.  )  (See  also  [5]  .  ) 


Proof.  In  any  partition  of  an  n-gon,  there  are  n-2  non-intersecting 
triangles  formed  by  n-3  arcs  and  n  sides.  And  for  any  n  2  4,  no 
triangle  can  be  formed  by  3  sides.  Let  x  be  the  number  of  triangles 
with  two  sides  and  one  arc,  y  be  the  number  of  triangles  with  one  side 
and  two  arcs,  and  7.  be  the  number  of  triangles  with  three  ares. 

Since  an  arc  is  used  in  two  triangles,  we  have 


x  +  2y  +  3z  =  2(n-3)  , 
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Since  the  polygon  has  n  sides,  we  have 

2x  +  y  =  n  .  (5) 

From  (4)  and  (5),  we  get 

3x  =  3z  +  6 

Since  z  s  0,  we  have  x  S  2  .  | 


Lemma  4.  Let  P  and  P'  both  be  n-gons  where  the  corresponding 
weights  of  the  vertices  satisfy  w.  <  w/  ,  then  the  cost  of  an  optimum  parti 
tion  of  p  is  less  than  or  equal  to  the  cost  of  an  optimum  partition  of  P'  . 


Proof.  Omitted.  ■ 


If  we  use  .... 

partitioning  the  k-gon  with  weights 
stated  as 


w.  )  to  mean  the  minimum  cost  of 
k 

w.  optimally,  Lemma  4  can  be 

i 


C(w  .  w2.  .  .  .  ,wk)  <  CfWj'.w^ . w')  if  w.  *w'  . 


We  say  that  two  vertices  arc  connected  in  an  optimum  partition 
if  the  two  vertices  are  connected  by  an  arc  or  if  the  two  vertices  are 
adjacent  to  the  same  side. 


a 


In  the  rest  of  the  paper,  we  shall  use  V  ,  V_ . V  to  denote 

\  c.  n 

vertices  which  are  ordered  according  to  their  weights,  i.e.  w  <  w_  <  •  •  •  £  w  . 

12  n 

To  facilitate  the  presentation,  we  introduce  a  tie-breaking  rule  for  vertices 
of  equal  weights. 

If  there  are  two  or  more  vertices  with  weights  equal  to  the  smallest 

weight  Wj  ,  we  can  arbitrarily  choose  one  of  these  vertices  to  be  the  vertex 

.  Once  the  vertex  is  chosen,  further  ties  in  equal  weights  are  resolved 

by  regarding  the  vertex  which  is  closer  to  in  the  clockwise  direction  to  be 

of  less  weight.  With  this  tie-breaking  rule,  we  can  unambiguously  label  the 

vertices  V, ,  V_,  .  .  .  ,  V  for  each  choice  of  V,.  A  vertex  V.  is  said  to  be 
1  2  n  1  i 

smaller  than  another  vertex  V.  ,  denoted  by  V  <  V  ,  either  if  w  <  w.  or  if 
-  j  x  j  i  J 

w^  =  Wj  and  i  <  j.  We  say  that  V.  is  the  smallest  vertex  in  a  subpolygon 

if  it  is  smaller  than  any  other  vertices  in  the  subpolygon. 

After  the  vertices  are  labeled,  we  define  an  arc  V.-V.  to  be  less  than 

1  J  - 

another  arc  V  -V 

P  q 

if  min(i,j)  <  min(p,  q) 

f  min(i,  j)  =  min(p,  q) 
or 

(  max(i,  j)  <  max(p,  q)  . 

(For  example,  the  arc  *s  ^ess  than  the  arc  V^-V^  .  )  Every  partition  of 

an  n-gon  has  n-3  arcs  which  can  be  sorted  from  the  smallest  to  the  largest 
into  an  ordered  sequence  of  arcs,  i.e.,  each  partition  is  associated  with  a 
unique  ordered  sequence  of  arcs.  We  define  a  partition  P  to  be  lexicographi- 
cally  less  than  a  partition  Q  if  the  ordered  sequence  of  arcs  associated  with 
P  is  lexicographically  less  than  that  associated  with  Q. 


When  there  is  more  than  one  optimum  partition,  we  use  the 
I  -optimum  partition  (i.  e.  ,  lexicographically-optimum  partition)  to  mean 
the  lexicographically  smallest  optimum  partition,  and  use  an  optimum  parti¬ 
tion  to  mean  some  partition  of  minimum  cost. 

We  shall  use  V  ,  V,  ,  . .  .  to  denote  vertices  which  are  unordered 
a  b 

in  weights,  and  T  to  denote  the  product  of  the  weights  of  any  three 
ijk 

vertices  V.,  V.  and  V,  . 

i  J  k 

Theorem  1.  For  every  way  of  choosing  V^.V^,  ...  (as  prescribed),  there 
is  always  an  optimum  partition  containing  V ^-V ^  and  V^-V^.  (Here,  V^-V^ 
and  V  -V  may  be  either  arcs  or  sides.  ) 

Proof:  The  proof  is  by  induction.  For  the  optimum  partitions  of  a  triangle 
and  a  4-go n,  the  theorem  is  true.  Assume  that  the  theorem  is  true  for  all 
k-gons  (3ikin-l)  and  consider  the  optimum  partitions  of  an  n-gon. 

From  Lemma  3,  in  any  optimum  partition,  we  can  find  at  least  two 
vertices  having  degree  two.  Call  these  two  vertices  W  and  V.  .  We  can 
divide  this  into  two  cases. 

(i)  One  of  the  two  vertices  V.  (or  V.)  is  not  V,,  V  or  V  in  some  optimum 

l  j  1  £  -i 

partition  of  the  n-gon.  In  this  case,  we  can  remove  the  vertex  V  with 
its  two  sides  and  obtain  an  (n-l)-gon.  In  this  (n-l)-gon,  V^,  V 
are  the  three  vertices  with  smallest  weights.  By  the  induction  assump¬ 
tion,  V  is  connected  to  both  V 2  and  V3  in  an  optimum  partition. 


(ii)  Consider  the  complementary  case  of  (i),  in  all  the  optimum  partitions  of 

the  n-gon,  all  the  vertices  with  degree  two  are  from  the  set  [V^.V^.  V^). 
(In  this  case,  there  will  be  at  most  three  vertices  with  degree  two  in  every 
optimum  partition,  )  We  have  the  following  three  subcases: 

(a)  V.  =  V_  and  V.  =  V  in  some  optimum  partition  of  the  n-gon, 

i  2  j  i 

i.e.  ,  both  and  have  degree  two  simultaneously.  In  this 
case,  we  first  remove  V 2  with  its  two  sides  and  form  an  (n-l)-gon. 

By  the  induction  assumption,  V  ,  must  be  connected  in  some  opti¬ 
mum  partition.  If  V  -V^  appears  as  an  arc,  it  reduces  to  (i).  So 

V  -  V  must  appear  as  a  side  of  the  (n-l)-gon,  and  reattaching  V  to 
the  (n-l)-gon  shows  that  either  V^,  and  are  mutually  adjacent 
or  V  -V3  is  a  side  of  the  n-gon.  In  the  former  case,  the  proof  is 
complete,  so  we  assume  that  V1'V3  is  a  side  the  n-gon.  Simi¬ 
larly,  we  can  remove  with  its  two  sides  and  show  that  V^, 

are  connected  by  a  side  of  the  n-gon. 

(b)  V  =  V  and  V  =  V_  in  some  optimum  partition  of  the  n-gon, 

i  1  J  2 

i.e.  ,  V  and  V_  both  have  degree  two  simultaneously.  In  this 

1  2 

case,  we  can  first  remove  and  form  an  (n-l)-gon  where  V^»  V^, 

V  are  the  three  vertices  with  smallest  weights.  By  the  induction 
4 

assumption,  V ^  is  connected  to  both  and  in  an  optimum 

partition.  If  V^-V^  or  ^2*^4  aPPears  as  an  arc>  ii  reduces  to  (i;. 

Hence,  V_-V  and  V_-V  must  both  be  sides  of  the  n-gon.  Simi- 

larly,  we  can  remove  V 2  with  its  two  sides  and  form  an  (n-l)-gon 

where  V  ,  V  ,  V  are  the  three  vertices  with  smallest  weights. 

13  4 
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Again,  must  be  connected  to  and  by  6ides  of  the  n-gon. 
But  for  any  n-gon  with  n  ^  5,  it  is  impossible  to  have  and 
both  adjacent  to  and  V ^  at  the  same  time,  i.e.  ,  and 
cannot  both  have  degree  two  in  an  optimum  partition  of  any  n-gon 
with  n  ^  5. 

(c)  V.  =  V^,  V.  =  in  some  optimum  partition  of  the  n-gon.  By 
argument  similar  to  (b),  we  can  show  that  must  be  adjacent 
to  V  and  in  the  n-gon.  The  situation  is  as  shown  in  Fig.  3(a) 
Then  the  partition  in  Fig.  3(b)  is  cheaper  because 


T  ^  T 
123  I2q 


and  C(w1,wq,wy,wt,wx,vw3)<  C(w2’  Wq’  Wy’  Wt’  Wx’  Wp’  W3) 


according  to  Lemma  4. 


Corollary  1.  For  every  way  of  choosing  V^.V^,  ...  (as  prescribed),  the 
t  -optimum  partition  always  contains  V  -V  and  V  -V  . 

\  Cm  1  J 


Proof:  It  follows  from  Theorem  1  and  the  definition  of  the  t  -optimum 
partition.  B 


Once  we  know  V  -V.  and  V  -V  always  exist  in  the  I  -optimum  parti- 

1  M  1  J 

tion,  we  can  use  this  fact  recursively.  Hence,  in  finding  the  /  -optimum 
partition  of  a  given  polygon,  we  can  decompose  it  into  subpolygons  by  joining 
the  smallest  vertex  with  the  second  smallest  and  third  smallest  vertices 
repeatedly,  until  each  of  these  subpolygcns  has  the  property  that  its  smallest 
vertex  is  adjacent  to  both  its  second  smallest  and  the  third  smallest  vertices. 

A  polygon  having  adjacent  to  and  by  sides  will  be  called 
a  basic  polygon. 


Theorem  2.  A  necessary  but  not  sufficient  condition  for  V^-V^  to  exist  in  an 
optimum  partition  of  a  basic  polygon  is 


_L  +  J_  s  ±  +JL 

w,  w .  w,  W- 
14  2  3 


Furthermore,  if  V ^-V ^  is  not  present  in  the  f  -optimum  partition, 
then  V  ,  V4  are  always  connected  in  the  /  -optimum  partition. 

Proof.  If  V2’ V3  are  not  connected  in  the  1  -optimum  partition  of  a  basic 
polygon,  the  degree  of  is  greater  than  or  equal  to  3.  Let  be  a  vertex 
in  the  polygon  and  vp  are  connected  in  the  f  -optimum  partition.  V ^  is 
either  in  the  subpolygon  containing  V.,  V  and  V  or  in  the  subpolygon  con- 

i  c  p 

taining  V  ,  V,  and  V  .  In  either  case,  V  will  be  the  third  smallest  vertex 
1  3  p  4 

in  the  subpolygon.  From  Corollary  1,  V  ,  V4  are  connected  in  the  I -optimum 
partition  of  the  subpolygon  and  it  also  follows  that  V  ,  are  connected  in  the 
/-optimum  partition  of  the  basic  polygon. 
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If  V  ,  V  are  connected  in  an  optimum  partition,  then  we  have  an 
2  3 

(n-l)-gon  where  V  is  the  smallest  vertex  and  V4  is  the  third  smallest 

vertex.  By  Theorem  1,  there  exists  an  optimum  partition  of  the  (n-l)-gon 

in  which  V,  ,  V.  are  connected.  Thus  by  induction  on  n,  we  can  assume 
2  4 

that  V ^  is  adjacent  to  in  the  basic  polygon  as  shown  in  Fig.  4(a). 


Fig.  4 


The  cost  of  the  partition  in  Fig.  4(a)  is 


Tj  23  +  C(w2»  w^,  . .  .  ,  w^t  . . .  >  w^)  . 

And  the  cost  of  the  partition  in  Fig.  4(b)  is 

Tl24  +  C(Wj<w4 . wt<....w3)  . 


(7) 


(8) 
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According  to  Lemma  4, 


c(wrw4,  ...,w3)  <  C(w2,w4 . wt,...,w3)  .  (9) 

Since  the  weights  of  the  vertices  between  V4  and  in  the  clockwise  direction 
are  all  greater  than  or  equal  to  w4>  the  difference  between  RHS  and  LHS 


of  (9)  is  at  least 


T  -  T 
243  143  * 


So  the  necessary  condition  for  (7)  to  be  no  greater  than  (8)  is 


or 


T123+T243  "T124+T134 


_L+  _L  £  _L  +  _L 

W1  w4  w2  W3 


Lemma  5.  In  an  optimum  partition  of  an  n-gon,  let  V  V  ,V  and  V  be 

x’  y*  z  w 

four  vertices  of  an  inscribed  quadrilateral  (V^  and  are  not  adjacent  in 
the  quadrilateral).  A  necessary  condition  for  V  -V  to  exist  is 


— L+  J_  __L  +  1 

w  w  w  w 
x  z  y  w 


(10) 


Proof:  The  cost  of  partitioning  the  quadrilateral  by  the  arc  V  -V  is 

x  z 


T  +  T 
xyz  xzw 


(ID 


and  the  cost  of  partitioning  the  quadrilateral  by  the  arc  V  -V  is 

y  w 

T  +  T 
xyw  yzw 


(12) 


For  optimality,  we  have  (11)  £  (12)  which  is  (10). 
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Note  that  if  strict  inequality  holds  in  (10),  the  necessary  condition  is 

also  sufficient.  If  equality  holds  in  (10),  the  condition  is  sufficient  for  V  -V 

X  z 

to  exist  in  the  1  -optimum  partition  provided  min(x,  z)  <  min(y,w).  This  lemma 

is  a  generalization  of  Lemma  1  of  Chin  [3]  where  is  the  vertex  with  the 

smallest  weight  and  V  ,  V  ,  V  are  three  consecutive  vertices  with  w 

greater  than  both  w  and  w  . 

6  x  z 

A  partition  is  called  stable  if  every  quadrilateral  in  the  partition 
satisfies  (10). 


Corollary  2.  An  optimum  partition  is  stable  but  a  stable  partition  may  not 
be  optimum. 


Proof.  The  fact  that  optimum  partition  has  to  be  stable  follows  from  Lemma  5. 
Figure  5  gives  an  example  that  a  stable  partition  may  not  be  optimum.  ■ 


Fig-  5 


In  any  partition  of  an  n-gonf  every  arc  dissects  a  unique  quadri¬ 
lateral  Let  V  V  V  V  be  the  four  vertices  of  an  inscribed  quadri- 
'  x '  y '  z  w 

lateral  and  V  -V  be  the  arc  which  dissects  the  quadrilateral.  We  define 
X  z 

V  -V  to  be  a  vertical  arc  if  (13)  or  (14)  is  satisfied, 
x  z  - 
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min(w  ,  w  )  <  ntin(w  ,  \v  1 
X  z  y  w 


min(w  ,  w  )  =  min(w  ,  w  1 

x  z  y  w  | 

max(w  ,  w  )  <  max(w  ,  \v  )  I 
x  z  y  w  ) 


We  define  V  -V  to  be  a  horizontal  are  if  (15)  is  satisfied 
X  z - 


min(w  ,  w  )  >  min(w  ,  w  ) 
X  z  y  w 


max(w  ,w  )<  max(w  ,w  ) 
X  Z  y  W 


For  brevity,  we  shall  use  h-arcs  and  v-arcs  to  denote  horizontal  arcs 


and  vertical  arcs  from  now  on. 


Corollary  3.  All  arcs  in  an  optimum  partition  must  be  cither  vertical 


arcs  or  horizontal  arcs. 


Proof:  Let  V  -V  be  an  arc  which  is  neither  vertical  nor  horizontal. 
-  x  z 

There  are  two  cases: 

Case  1.  min(w  ,  w  )  =  min(w  ,  w  ) 

x  z  y  w 

and  max(w  ,  w  )  >  max(w  ,  w  ) 

x  z  y  w 

Case  2.  minfw  ,  w  )  >  mi n(w  ,  w  ) 

x  7,  y  w 

and  max(w  .  w  )  >-  mnx(w  ,  w  )  . 

x  z  y  w 

In  both  cases,  the  inequality  (10)  in  Lemma  5  cannot  be  satisfied. 


This  implies  that  the  partition  is  not  stable  and  hence  cannot  ho  optimum. 
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Theorem  3.  Fet  and  be  two  arbitrary  vertices  winch  arc  not  adjacent 

in  a  polygon,  and  V  be  the  smallest  vertex  from  V  to  V  in  the  clockwise 
w  x  z 

manner  (V  ^  V  V  /  V  )  and  V  be  the  smallest  vertex  from  V  to  V 
w  x  w  z  y  z  x 

in  the  clockwise  manner  (V  ^  V  ,  V  f  V  ),  This  is  shown  in  Fig.  6  where 

Y  x  y  7, 


without  loss  of  generality,  we  assume  that  V  <  V  and  V  <  V  .  A 

x  z  y  w  ’ 

necessary  condition  for  V  -  V  to  exist  as  an  h-arc  in  the  I  -optimum 

X  z 

partition  is  that 


w  <  w  <  w  <  w 
y  x  z  w 


(Note  that  the  necessary  condition  still  holds  when  the  positions  of  V  and  V 

y  w 


are  interchanged.  ) 


.  •-  Vw 

„  — {Ww  V  _ 

#■  \  -  > 

/ 

/  \ 

Vx'Wx) 


Vy 


Fig.  6 


Proof.  The  proof  is  by  contradiction.  If  w  <:  w  ,  w  must  be  equal  to 

x  y  x  1 

the  smallest  weight  w  and  V  -V  can  never  satisfy  (15).  Hence,  in 

1  X  z 

order  that  vx‘v2  exists  as  an  h-arc  in  the  f  -optimum  partition,  we  must 
have  w  <w  <  w  .  Since  V  is  the  smallest  vertex  from  V  to  V  in 

y  x  y  z  x 

the  clockwise  manner  and  V  <  V  ,  we  must  have  V  =  V 

x  w  y  1  ’ 
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Assume  for  the  moment  that  V  ,  <  V  <  V  .  From  Corollary  1, 

3  x  z 

both  V  -V  and  V  -V  exist  in  the  i  -optimum  partition,  and  the  two  arcs 

\  C,  1  J 

would  divide  the  polygon  into  subpolygons.  If  and  are  in  different 

subpolygons,  then  they  cannot  be  connected  in  the  l  -optimum  partition. 

Without  loss  of  generality,  we  can  assume  that  the  polygon  is  a  basic  polygon. 
In  this  basic  polygon,  either  or  ^i"^4  ex^s^s  the  ^  -optimum 

partition  (Theorem  2). 

If  V.,  V,  are  connected,  then  V  and  V  are  both  in  a  smaller  polygon 
C  j  X  z  * 

in  which  we  can  treat  V  as  the  smallest  vertex  and  repeat  the  argument.  If 

V  ,  are  connected,  the  basic  polygon  is  again  divided  into  two  subpolygons 

and  V  and  V  both  have  to  be  in  one  of  the  subpolygons  and  the  subpolygon  has 
X  z 

at  most  n-1  sides.  (Otherwise  V  -V  can  never  exist  in  the  t  -optimum 

X  z 

partition.  )  The  successive  reduction  in  the  size  of  the  polygon  will  either 

make  the  connection  V  -V  impossible,  or  force  V  and  V  to  become  the 

x  z  x  z 

second  smallest  and  the  third  smallest  vertices  in  a  basic  subpolygon.  Let  V 

m 

be  the  smallest  vertex  in  this  basic  subpolygon.  In  order  that  V  -V  appear 

X  z 

as  an  h-arc,  we  must  have  w  >  w  .  From  Theorem  2,  the  necessary  condi- 

x  m 

tion  for  (i.  e.  ^2*^3^  to  ex*st  *n  an  optimum  partition  of  the  subpolygon 

is 

-!+  JL  a  +  . 

w  w  w  w 
x  z  m  w 


Since  w  >  w  ,  the  inequality  is  valid  only  if  w  <  w 
x  m  1  ■* 


Corollary  4.  A  weaker  necessary  condition  for  to  exist  as  an  h-arc 


in  the  t  -optimum  partition  is  that 


V  <  V  <  V  <  V 
y  x  z  w 


Proof.  This  follows  from  Theorem  3.  ■ 
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We  call  any  arc  which  satisfies  this  weaker  necessary  condition  a 
potential  h-arc.  Let  P  be  the  set  of  potential  h-arcs  in  the  n-gon  and  H 
be  the  set  of  h-arcs  in  the  i -optimum  partition,  we  have  P  3  H  where  the 
inclusion  could  be  proper. 

Corollary  5.  Let  be  the  largest  vertex  in  the  polygon  and  and 
be  its  two  neighboring  vertices.  If  there  exists  a  vertex  V  such  that 

y 

V  <  V  and  V  <  V  ,  then  V  -V  is  a  potential  h-arc. 
y  x  y  z  x  z  r 

Proof.  This  follows  directly  from  Corollary  4  where  there  is  only  one  vertex 

between  V  and  V  .  H 
x  z 

Two  arcs  are  called  compatible  if  both  arcs  can  exist  simultaneously 

in  a  partition.  Assume  that  all  weights  of  the  vertices  are  distinct,  then  there 

are  (n-1)!  distinct  permutations  of  the  weights  around  an  n-gon.  For 

example,  the  weights  10,  11,  25,  40,  1 2  in  Fig.  5(a)  correspond  to  the 

permutation  w  ,w,,w  ,w  ,w,  (where  w  <  w,  <  w_  <  w  <  w_  ).  There 
\  4  4  5  3  1  c  3  4  5 

are  infinitely  many  values  of  the  weights  which  correspond  to  the  same  per¬ 
mutation.  For  example,  1,  16,  34,77,  29  also  corresponds  to  w, ,  w_,  w.,  w_,  w 

1  c  4  D 

but  its  optimum  partition  is  different  from  that  of  10,  11,  25,  40,  12.  However, 
all  the  potential  h-arcs  in  all  the  n-gons  with  the  same  permutation  of  weights 
are  compatible.  We  state  this  remarkable  fact  as  Theorem  4. 

Theorem  4.  All  potential  h-arcs  arc  compatible. 

Proof.  The  proof  is  by  contradiction.  Let  V  ,  V  ,  V  and  V  be  the  four 
-  x  y  z  w 

vertices  described  in  Theorem  3.  Hence,  we  have  V  <  V  <  V  <  V 

y  x  z  w 


and  V  -V  is  a  potential  h-arc.  Let  V  -V  be  a  potential  h-arc  which  is  not 
x  Z  PH 

compatible  to  V  -V  ,  as  shown  in  Fig.  7.  Without  loss  of  generality,  we 
r  x  z 

can  assume  V  <  V  .  (The  proof  for  the  case  V  <  V  is  similar  to  that 


P  q 


which  follows.  ) 


q  P 


%  o 


I 


-  - 


Fig.  7 


Since  V  is  the  smallest  vertex  between  V  and  V  in  the  clockwise 
w  x  z 

manner,  we  have  V  <  V  <  V  .  Hence,  we  have  either  V  <  V  <  V  <  V 

z  w  q  y  p  z  q 

orV  <V  <V  <  V  .  Both  cases  violate  Corollary  4  and  V  -  V  cannot 
y  z  p  q  P  q 

be  a  potential  h-arc.  ■ 

Note  that  the  potential  h-arc  V  -V  always  dissects  the  n-gon  into  two 

X  z 

subpolygons  and  one  of  these  subpolygons  has  the  property  that  all  its  vertices 
except  V  and  V  have  weights  no  smaller  than  max(w  ,w  ).  Wc  shall  call  this 

r  x  Z  X  z 

subpolygon  the  upper  subpolygon  of  vx_vz  •  For  example,  the  subpolygon 

V  . V  . V  . V  in  Fig.  7  is  the  upper  subpolygon  of  V  -V  . 

x  w  q  z  x  z 
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Using  Corollary  4  and  Theorem  4,  we  can  generate  all  the  potential 


h-arcs  of  a  polygon. 

Let  V  -V  be  the  arc  defined  in  Corollary  5,  i.e.  V,  <  V  <  V  <  V 
x  z  1  x  z  w 

The  arc  V  -V  is  a  potential  h-arc  compatible  to  all  other  potential  h-arcs  in 
X  z 

the  n-gon.  Furthermore,  there  is  no  other  potential  h-arc  in  its  upper  subpoly 

gon.  Now  consider  the  (n-l)-gon  obtained  by  cutting  out  V  .  In  this  (n-1)- 

gon,  let  V  ,be  the  largest  vertex  and  V  ,  and  V  ,  be  the  two  neighbors  of 
w  x  z 

V  /  where  V  <  V  ,<  V  .<  V  /  .  Then  V  /-V  /  is  again  a  potential  h-arc 

W  1  X  Z  W  X  z 

compatible  to  all  other  potential  h-arcs  in  the  n-gon  and  there  is  no  other 
potential  h-arc  in  its  upper  subpolygon  which  has  not  been  generated.  This 
is  true  even  if  V  is  in  the  upper  subpolygon  of  ,  -V^  , .  If  we  repeat  the 
process  of  cutting  out  the  largest  vertex,  we  get  a  set  P  of  arcs,  all  arcs 
satisfy  Corollary  4.  The  h-arcs  of  the  i -optimum  partition  must  be  a 
subset  of  these  arcs. 

The  process  of  cutting  out  the  largest  vertex  can  be  made  into  an 
algorithm  which  is  O(n).  We  shall  call  this  algorithm  the  one -sweep 
algorithm.  The  output  of  the  one- sweep  algorithm  is  a  set  S  of  n-3  arcs. 

S  is  empty  initially. 

The  one- sweep  algorithm: 

Starting  from  the  smallest  vertex,  say  V^,  we  travel  in  the  clockwise  direc¬ 
tion  around  the  polygon  and  push  the  weights  of  the  vertices  successively  onto 
the  stack  as  follows  (w^  will  be  at  the  bottom  of  the  stack). 

(a)  Let  Vt  be  the  top  element  on  the  stack,  be  the  element  immedi¬ 

ately  below  V  ,  and  Vc  be  the  element  to  be  pushed  onto  the  stack. 


If  there  are  two  or  more  vertices  on  the  stack  and  w  >w  ,  add 

t  c 

V,  ,  -V  to  S,  pop  V,  off  the  stack;  if  there  is  only  one  vertex  on 
t-1  c  r  t 

the  stack  or  w^  £  wc  ,  push  wc  onto  the  stack.  Repeat  this  step 
th 

until  the  n  vertex  has  been  pushed  onto  the  stack. 

(b)  If  there  are  more  than  three  vertices  on  the  stack,  add  V  -V 

t-1  c 

to  S,  pop  V  off  the  stack  and  repeat  this  step,  else  stop. 

Since  we  do  not  check  for  the  existence  of  a  smallest  vertex  whose 
weight  is  strictly  no  larger  than  those  of  the  two  neighbors  of  the  largest  ver 
tex,  i.  e.  the  existence  of  the  vertex  V  in  Corollary  4,  not  all  the  n-3  arcs  g 

y 

erated  by  the  algorithm  are  potential  h-arcs.  However,  it  is  not  difficult  to 
verify  that  the  one-sweep  algorithm  always  generates  a  set  S  of  n-  3  arcs 
which  contains  the  set  P  of  all  potential  h-arcs  which  contains  the  set  H  of 
all  h-arcs  in  the  t  -optimum  partition  of  the  n-gon,  i.  e.  , 

S  2  P  2  H 

where  each  inclusion  could  be  proper.  For  example,  if  the  weights  of  the 

vertices  around  the  n-gon  in  the  clockwise  direction  are  w  ,  w0 . w 

1  c  n 

where  w,  ^  w.,  <  •  •  •  ^  w  ,  none  of  the  arcs  in  the  n-gon  can  satisfy 
1  2  n  ’ 

Corollary  4  and  hence  there  are  no  potential  h-arcs  in  the  n-gon.  The  one- 
sweep  algorithm  would  still  generate  n-3  arcs  for  the  n-gon  but  none  of  the 
arcs  generated  is  a  potential  h-arc. 


3.  Conclusion 


In  this  paper,  we  have  presented  several  theorems  on  the  Polygon 
Partitioning  Problem.  Some  of  these  theorems  are  characterizations  of  the 
optimum  partitions  of  any  n- sided  convex  polygon,  while  the  others  apply 
to  the  unique  lexicographically  smallest  optimum  partition.  Based  on  these 
theorems,  an  O(n)  algorithm  for  finding  a  near- optimum  partition  can  be 
developed  1 12] .  The  cost  of  the  partition  produced  by  the  heuristic  algorithm 
never  exceeds  1.  155  Copt,  where  Copt  is  the  optimum  cost  of  partitioning 
the  polygon.  An  0(n  log  n)  algorithm  for  finding  the  unique  lexicographically 
smallest  optimum  partition  will  be  presented  in  part  II. 
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Abstract: 


This  paper  considers  the  computation  of  matrix  chain  products 

of  the  form  M.  X  X  ...  X  M  ,  .  If  the  matrices  arc  of  different 
12  n- 1 

dimensions,  the  order  in  which  the  matrices  are  computed  affects  the 
number  of  operations.  An  optimum  order  is  an  order  which  minimizes 
the  total  number  of  operations.  Some  theorems  about  an  optimum  order 
of  computing  the  matrices  have  been  presented  in  part  I.  Based  on 
those  theorems,  an  0(n  log  n)  algorithm  for  finding  the  optimum  order 
is  presented  here. 
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1.  Introduction 


In  Part  I  of  this  paper,  we  have  transformed  the  matrix  chain 
product  problem  into  the  optimum  partitioning  problem  and  have  stated 
several  theorems  about  the  optimum  partitions  of  an  n- sided  convex  polygon. 
Based  on  these  theorems,  we  now  present  algorithms  for  finding  the  unique 
(  -optimum  (lexicographically  smallest  optimum)  partition. 

Using  the  same  notation  as  in  Part  I  of  this  paper,  we  can  assume 

that  we  have  uniquely  labelled  all  vertices  of  the  n-gon.  A  partition  is 

called  a  fan  if  it  consists  of  only  v-arcs  joining  the  smallest  vertex  to  all 

other  vertices  in  the  polygon.  We  shall  denote  the  fan  of  a  polygon 

V  -V,  -V  ....  -V  by  Fan(  w  I  w,  .  w  ,  .  .  .  ,  w  ).  The  smallest  vertex  V 
lbc  n  r  1  b  c  n  1 

is  called  the  cente r  of  the  fan. 

We  define  a  vertex  as  a  local  maximum  vertex  if  it  is  larger  than  its 
two  neighbors  and  define  a  vertex  as  a  local  minimum  vertex  if  it  is  smaller 
than  its  two  neighbors.  A  polygon  is  called  a  monotone  polygon  if  there 
exists  only  one  local  maximum  and  one  local  minimum  vertex.  We  shall 
first  give  an  O(n)  algorithm  for  finding  the  l  -optimum  partition  of  a  mono¬ 
tone  polygon  and  then  give  an  0(n  log  n)  algorithm  for  finding  the  i  -optimum 
partition  of  a  general  convex  polygon. 


2.  Monotone  Basic  Polygon 

In  this  section,  let  us  consider  the  optimum  partition  of  a  monotone 
polygon,  i.e,  a  polygon  with  only  one  local  minimum  vertex  and  one  local 
maximum  vertex.  It  follows  from  Corollary  1  of  Part  I  that  we  can 
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consider  a  monotone  basic  polygon  only.  The  understanding  of  this  special  case 


) 


I 

l 


is  necessary  in  finding  the  optimum  partition  of  a  general  convex  polygon. 
Consider  a  monotone  basic  n-gon  V  -V  -V  -  •  •  •  -V  ,  the  fan  of  the 

1  C*  C  j 

polygon  is  denoted  by 

Fan(w.  wvw  ,  ...,w,) 

1  c.  c  5 

where  the  smallest  vertex  is  the  center  of  the  fan. 

The  definition  of  a  fan  can  also  be  applied  to  subpolygons  as  well.  For 
example,  if  V^,  are  connecte(i  in  the  basic  n-gon  and  becomes  the 
smallest  vertex  in  the  (n-l)-sided  subpolygon,  the  partition  formed  by  con¬ 
necting  to  all  vertices  in  the  (n-l)-gon  is  denoted  by 

Fan(w^|  w^P  . .  .  ,  w^)  . 

Lemma  1.  If  none  of  the  potential  h-arcs  appears  in  the  £  -optimum  parti¬ 
tion  of  the  n-gon,  the  £  -optimum  partition  must  be  the  fan  of  the  n-gon. 

Proof.  From  Theorem  3  of  Part  I,  we  know  that  any  arc  which  exists  as 
an  h-arc  in  the  (  -optimum  partition  mus  je  a  potential  h-arc.  Hence,  if 
the  i  -optimum  partition  does  not  contain  any  potential  h-arc,  the  £  -optimum 
partition  must  be  made  up  of  v-arcs  only.  Hence,  wc  lave  to  show  that 
among  all  partitions  which  are  made  up  of  v-arcs  only,  the  fan  is  (i)  the 
lexicographically  smallest  and  (ii)  one  of  the  cheapest  partitions  in  the  n-gon. 

(i)  Since  the  fan  consists  of  only  v-arcs  joining  to  all  other  vertices 
in  the  n-gon,  it  is  by  definition  the  lexicographically  smallest  partition. 

(ii)  Suppose  the  t  -optimum  partition  contains  v-arcs  only  hut  is  not  the 
fan.  There  must  exist  three  vertices  V^,  V^,  V.  such  that  the  triangles 
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V  V  V.  and  V  V  V  are  present  in  the  l  -optimum  partition.  Since  V.-V. 
1  1  j  i  j  k  1  j 

is  a  v-arc  (by  assumption)  and  is  the  smallest  vertex  in  the  n-gon,  we 


have  w,  =  min(w.,w.)  and  max(w. ,  w. )  ^  w,  .  If  we  replace  the  v-arc 
1  i  j  i  J  k 

V  -V  by  the  v-arc  V  -V  ,  we  can  get  a  partition  whose  cost  is  less  than  or 
i  j  lk 

equal  to  that  of  the  (  -optimum  partition  but  is  lexcographically  smaller 


than  the  t -optimum  partition,  and  results  in  a  contradiction.  ■ 


Let  V.-V.  and  V  -V  be  two  potential  h-arcs  of  any  n-gon.  We 

i  j  p  q 

say  that  V.-V.  is  above  V  -V  (and  V  -V  is  below  V.-V.)  if  the  upper 
l  j  -  p  q  p  q  -  l  j 

subpolygon  of  V^-V  contains  the  upper  subpolygon  of  V.-V.  . 

Let  P  be  the  set  of  all  potential  h-arcs  in  a  monotone  basic  n-gon. 


P  can  have  at  most  (n-3)  arcs. 


Lemma  2.  For  any  two  arcs  in  P,  say  V.-V.  and  V  -V  ,  we  must  have 
-  y  x  •,  p  q 

either  V.-V.  above  V  -V  or  V  -V  above  V.-V.  . 

1  J  P  q  P  q  *  J 

Proof.  By  contradiction.  Let  V.-V.  and  V  -V  be  two  arcs  in  P  which 

-  i  j  p  q 

do  not  satisfy  this  lemma.  Then  the  intersection  of  the  upper  subpolygons 

of  V.-V.  and  V  -V  must  either  be  empty  or  consists  of  part  of  each  upper 
i  j  p  q 

subpolygon  only. 

Since  the  vertices  other  than  V.,  V  in  the  upper  subpolygon  of 

V.-V.  must  have  weights  larger  than  max(w. ,w.),  the  local  maximum  vertex 
1  J  6  6  1  J 

of  the  monotone  basic  polygon  must  be  present  in  the  upper  subpolygon  of 


V.-V  .  Similarly,  the  local  maximum  vertex  of  the  monotone  basic  polygon 
must  also  be  present  in  the  upper  subpolygon  of  V^-V^.  Hence,  the  inter¬ 
sections  of  the  upper  subpolygons  of  V.-V  and  V  -V^  cannot  be  empty. 
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From  Theorem  4  of  Part  I,  we  know  that  V.-V.  and  V  -V  cannot 

1  J  P  q 

cross  each  other  and  hence  the  intersection  of  their  upper  subpolygons 

cannot  consist  of  part  of  each  upper  subpolygons  only.  B 

We  can  actually  show  this  ordering  of  potential  h-arcs  pictorially 

by  drawing  a  monotone  basic  polygon  in  such  a  way  that  the  local  maximum 

vertex  is  always  at  the  top  and  the  local  minimum  vertex  is  at  the  bottom. 

Then  a  potential  h-arc  V.-V.  is  physically  above  another  potential  h-arc 

V  -V  if  the  upper  subpolygon  of  V  -V  contains  the  upper  subpolygon  of 
p  q  p  q 

V.-V.  .  From  the  definition  of  the  upper  subpolygon,  we  can  see  that 
min(w.,w.)>  max(w  ,  w  )  if  V.-V.  is  above  V  -V 

1  j  p  q  1  j  p  q 

Consider  the  monotone  basic  n-gon  which  is  shown  symbolically  in 

Figure  1.  V  is  the  local  maximum  vertex  and  V.-V.,  V  -V  are  potential 

«  1  J  P  q 

h-arcs  of  the  monotone  basic  n-gon.  The  subpolygon  V  -  •  •  •  -  V.-V. . V 

P  i  J  q 

which  is  formed  by  tv  o  potential  h-arcs  V  -V  and  V.-V.  and  the  sides  of 

p  q  1  J 

the  n-gon  from  V  to  V.  and  from  V.  to  V  in  the  clockwise  direction  is  said 

p  i  j  q 

to  be  bounded  above  by  the  potential  h-arc  V.-V.  and  bounded  below  by  the 


potential  h-arc  V  -V 

P  q 
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Figure  1 


Lemma  3.  Any  subpolygon  which  is  bounded  by  two  potential  h-arcs  of  the 
monotone  basic  n-gon  is  itself  a  monotone  polygon. 


Proof.  Consider  the  subpolygon  V  -  •  •  •  -V.-V •  •  •  -V  in  Figure  1. 

-  P  i  i  q 


Without  loss  of  generality,  we  can  assume  V.  <  V.  and  V  <  V  .  Since 

i  J  p  q 


V  is  the  only  local  maximum  vertex  in  the  monotone  basic  n-gon,  we  must 


have  V,  <  V  < . .  •  <  V.  <  V  and  V  >  V.  >  . . .  >  V  >  V,  .  Hence,  V 
1  p  in  n  j  q  l 


is  the  unique  local  minimum  vertex  and  V.  is  the  unique  local  maximum 


vertex  in  the  subpolygon  -  ...  -  V.-V.  -  .  . .  -  ,  By  definition, 


-V.-V.-.-.-V  is  a  monotone  polygon. 


Lemma  4.  Any  potential  h-arc  of  a  subpolygon  bounded  above  and  below 
by  two  potential  h-arcs  of  the  monotone  basic  n-gon  is  also  a  potential 
h-arc  of  the  monotone  basic  n-gon. 

Proof.  Consider  the  subpolygon  V  -  • .  •  -V.-V.-  •  •  •  -V  in  Figure  1.  Let 

-  P  i  J  q 

V  -V  be  a  potential  h-arc  in  this  subpolygon  and  V  is  the  smallest  vertex 
X  z  w 

between  and  V^_  in  the  clockwise  direction  around  the  subpolygon.  With¬ 
out  loas  of  generality,  we  can  assume  V.  <V.  ,  V  <  V  and  V  <  V 
'  l  j  p  q  x  z 

Since  V  is  in  the  upper  subpolygon  of  the  potential  h-arc  V  -V  ,  we  have 
x  r  P  q 

Wj  <  w^  £  <  w^  <  w  .  Since  V\  <  any  vortex  in  the  upper  subpolygon 

of  V.-V.  and  V  <  V.  <  V  V  is  the  smallest  vertex  between  V  and  V  in 
i  J  w  l  j  w  x  z 

clockwise  direction  around  the  monotone  basic  n-gon.  Hence,  we  have 

w,  <  w  <  w  <  w  and  V  -V  is  a  potential  h-arc  of  the  monotone  basic 
1  x  z  w  x  z 

n-gon.  ■ 

We  can  now  summarize  what  we  have  discussed.  If  there 

is  no  h-arc  in  the  t -optimum  partition  of  a  monotone  basic 

n-gon,  the  t  -optimum  partition  must  be  a  fan.  Otherwise,  the  h-arcs  in  tne 

l  -optimum  partition  are  all  layered,  one  above  another.  If  we  consider  the 

local  maximum  vertex  V  and  the  local  minimum  vertex  V,  as  two  degen- 

n  1  ° 

erated  h-arcs,  then  the  l  -optimum  partition  of  a  monotone  basic  n-gon  will 
contain  one  or  more  monotone  subpolygons,  each  bounded  above  and  below 
by  two  h-arcs  and  the  f  -optimum  partition  of  each  of  these  monotone  sub¬ 
polygons  is  a  fan. 
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Then,  in  finding  the  i  -optimum  partition  of  a  monotone  basic  polygon,  we 
have  only  to  consider  those  partitions  which  contain  one  or  more  subpoly¬ 
gons  bounded  above  and  below  by  potential  h-arcs  and  each  of  these  sub¬ 
polygons  is  partitioned  by  a  fan.  Since  there  are  at  most  (n-3)  non¬ 
degenerated  potential  h-arcs  in  a  monotone  basic  n-gon,  there  will  be  at 
n-  3 

most  2  such  partitions  and  we  can  divide  all  these  partitions  into 

(n-2)  classes  by  the  number  of  non-degenerated  potential  h-arcs  a  partition 

contains.  These  classes  are  denoted  by  H  ,  H  ,...,  H  _  where  the  sub- 

U  1  n-3 

script  indicates  the  number  of  non-degenerated  potential  h-arcs  in  each 
partition  of  that  class. 

There  is  no  potential  h-arc  in  the  partitions  in  the  class  . 

Hence  the  class  consists  of  only  one  partition,  namely  the  fan 

Fan(w1 1  w^i  . .  .  ,  w^). 

In  the  class  H^,  each  partition  has  one  non-degenerated  potential 

h-arc.  Once  the  potential  h-arc  is  known,  the  rest  of  the  arcs  must  all  be 

vertical  arcs  forming  two  fans,  one  in  each  subpolygon. 

Two  typical  partitions  in  of  a  monotone  basic  polygon  are  shown 

in  Fig.  2.  In  Fig.  2a,  there  is  one  non-degcncrated  potential  h-arc, 

V  -V  (V  <  V.).  The  upper  subpolygon  is  a  fan 

Fan(w  I  w,,  . . .  ,  w. ) 
c 1  d  l 

and  the  lower  subpolygon  is  a  fan 

Fan(w,  |w,,w  ,w.,w-)  . 

1  c.  c  l  3 
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Fig.  2.  Two  typical  partitions  in  H  of  a  monoton 
10-gon.  1 


In  Fig.  2b,  there  is  one  potential  h-arc,  V ^-V ^  ,  and  the  upper 


subpolygon  is  a  fan 


Fan(w  |  w  ,  .  .  .  ,  wj 
c  c  5 


and  the  lower  subpolygon  is  a  degenerated  fan,  a  triangle. 

Assume  that  V^-V ^  is  the  only  h-arc,  then  the  cost  is  (see  Fig.  2b) 


w  w_w,  +  w  (w  w  ,  +  w  ,w  +  w  w,  +w,w  +  w  w,  +  w,  w  +w  w  ) 
123  2  c  d  de  ef  fg  gh  hi  i  3; 


=  T123  +W2(wc:W3)  ’ 


(1) 


where  w^w^  is  the  shorthand  notation  of  the  sum  of  adjacent  products  from 


w  to  w.  in  the  clockwise  direction, 
c  3 


Note  that  the  cost  of  H  of  the  polygon  shown  in  Fig.  2  is 


Fan(w1 1  w.,,  •  •  .  .  w^) 

=  Wl(w2:w3) 

The  condition  of  (1)  to  be  less  than  (2)  is 

w2.(wc:w3) 

(w2:w3)  -  w2-  w3  <  W1 


(2) 


Similarly,  the  condition  for  the  partition  in  Fig.  2a  to  be  less 


than  Hq  is 


w  •  (w  w. ) 
c  d  i 

- - - - <  w,  . 

(w  :  w.)  -  w  •  w.  1 

Cl  Cl 


(3) 


We  say  that  a  partition  is  said  to  be  t  -optimal  among  the  partitions 

in  a  certain  class  (or  several  classes)  if  it  is  the  lexicographically  smallest 

partition  among  all  the  partitions  with  minimum  cost  in  that  class  (or 

several  classes).  Hence,  the  l  -optimum  partition  is  l  -optimal  among  all 

partitions  in  the  classes  H„,H . and  H  .  . 

0  l  n-3 

Now,  assume  that  the  f -optimal  partition  among  all  the  partitions  in 

H,,H, . H  ,  contains  only  one  potential  h-arc  V. -V,  only,  as  shown 

1  2  n-3  l  k 

in  Fig.  3.  (Note  that  V.-V^  will  exist  in  this  partition  as  an  h-arc.  )  This 
partition  will  be  the  l  -optimum  partition  of  the  monotone  basic  n-gon  if  it 
costs  less  than  that  of  the  fan  in  H^,  The  condition  that  the  partition  with 
V.-Vk  as  the  single  h-arc  costs  less  than  is 


w.  •  (w.  :  w,  ) 
l  i  k 

7 - -  <  w 

(w.  :  w,  )  -  w.  •  vv,  1 

l  k  i  k 


if  w.  w, 
i  k 


or 
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Fig.  3.  A  monotone  polygon  with  a  single  h-arc. 


where  C(w  ,  ,  .  ,  ,  w,  )  denotes  the  cost  of  the  optimum  partition  of  the 
1  k 

subpolygon  w. -w.  -  •  ••  -Wg-Wk  an<^  eclual  to  the  cost  of  the  fan  in  this 
case. 

An  h-arc  V  -V,  which  divides  a  polygon  into  two  subpolygons  is 
1  k 

called  a  positive  arc  with  respect  to  the  polygon  if  (4)  is  satisfied,  i.c.  ,  the 
partition  with  the  arc  as  the  only  h-arc  and  a  fan  in  each  of  the  two  subpoly¬ 
gons  costs  less  than  the  fan  in  the  same  polygon.  Otherwise,  it  is  called  a 
negative  arc  with  respect  to  the  polygon. 
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When  an  n-gon  is  divided  into  subpolygons,  an  h-arc  is  defined 
as  positive  in  a  subpolygon  if  the  cost  of  partition  of  the  subpolygon  with 
the  h-arc  as  the  only  h-arc  is  less  than  the  fan  in  the  subpolygon. 

Let  us  consider  a  partition  with  two  h-arcs  as  shown  in  Fig.  4, 
and  assume  that  this  partition  is  i -optimal  among  all  partitions  in  the 


Fig,  4.  A  monotone  8-gon  with  two  h-arcs. 

If  is  positive  with  respect  to  the  subpolygon 

VrVi-Vp-W  t^lcn  I"*16  condition  analogous  to  (4)  is 

C(w. ,  w  ,  w  ,  w,  ) 

.  .  t  P  q  k  _ _ 

l  (w- :  w,  )  -  [  (w  :w)-w  ■  w  ]  }  -  w  •  w,  <'W1 
l  k  p  q  p  q  i  k 

If  Is  positive  with  respect  to  the  whole  polygon 

.  -  Vk  f  then  the  condition  is 
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(5b) 


C(w  w  w  w  ,  w  ,w  ,w,) 
i  P  r  n  s’  q  k 

(w.:  w  )  -  w.  .  w!  <  W1 

i  k  l  k 


Note  that  (5b)  implies  (5a). 

The  condition  for  the  arc  V  -V  to  be  positive  with  respect  to 

PH 

the  subpolygon  V  -V  -V  -V  -V  -V  -V  is 
l  p  r  n  s  q  k 


C  (w  .  w  .  w  ,  w  .  w  ) 
p  r’  n  s  q 

— j - : - -  <  min(w. ,  w,  ) 

(w  :  w  )  -  w  •  w  i  k 


(6a) 


P  q 


p  q 


If  the  arc  V^-V^  is  positive  with  respect  to  the  whole  polygon 


-V  -  V  -V  -V  -V  -  V  -V  it  must  satisfy  (6b). 
1  l  p  r  n  s  q  k 


C  (w  .  w  ,  w  ,  w  ,  w  ) 
p  r’  n’  s’  q 

(w  :  w  )  -  w  •  w  1 

p  q  p  q 


(6b) 


Since  <  min(w.,wk),  condition  (6b)  implies  (6a). 

Here,  the  presence  of  V.-V  will  divide  the  original  polygon 

l  K 

into  two  subpolygons  where  V^- V  appears  in  the  upper  subpolygon. 

If  V^-V^  a  Pos’five  arc  with  respect  to  the  original  polygon,  then 
V  -V  is  certainly  positive  in  the  upper  subpolygon.  But  if  V  -V  is 

p  q  P  q 

positive  in  the  subpolygon,  the  arc  V  -V  may  become  negative  if 

p  q 

V.-V,  is  removed,  i.  e.  V  -V  becomes  negative  with  respect  to  the 
ik  P  q 

original  polygon. 


Similarly,  if  the  arc  V.-V^  is  positive  with  respect  to  a  sub- 
polygon,  the  arc  V^-V^  may  become  negative  if  the  arc  V_-V_  is 


p  q 


removed. 


.The  preceding  discussions  can  be  summarized  as  Theorem  1. 
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Theorem  1.  If  an  h-arc  is  positive  with  respect  to  a  polygon  then  the 

arc  is  positive  with  respect  to  any  subpolygon  containing  that  are.  If  an  h-arc 

is  positive  with  respect  to  a  subpolygon,  it  may  or  may  not  be  positive  with 

respect  to  a  larger  polygon  which  contains  the  subpolygon.  ■ 

There  are  two  intuitive  approaches  to  the  f  -optimum  partition  of 

a  monotone  basic  polygon.  The  first  approach  is  to  put  in  the  potential 

h-arcs  one  by  one.  Each  additional  potential  h-arc  will  improve  the  cost 

until  the  correct  number  of  h-arcs  is  reached.  Any  further  increase  in 

the  number  of  h-arcs  will  increase  the  cost.  To  introduce  an  h-arc  into 

the  polygon,  we  can  test  each  potential  h-arc  (at  most  n-3)  to  see  if  it 

is  positive  with  respect  to  the  whole  polygon.  If  yes,  that  positive 

arc  must  exist  in  the  l- optimum  partition,  and  the  polygon  will  be 

divided  into  two  subpolygons,  each  being  a  monotone  polygon.  We  can 

repeat  the  whole  process  of  testing  positiveness  of  the  h-arcs.  The 

trouble  is  that  all  these  arcs  may  be  negative  individually  with  respect  to 

the  whole  polygon  and  yet  may  not  be  the  optimum.  For  example,  two 

arcs  V.-V  an<i  may  negative  individually  with  respect  to  the 

whole  polygon  but  the  partition  with  both  V.-V.,  V  -V  present  at  the 

i  J  P  q 

same  time  may  cost  less  than  as  shown  in  Fig,  5a,  This  shows  that 
we  cannot  guarantee  an  optimum  partition  simply  because  no  more 
potential  h-arcs  can  be  added  one  at  a  time. 

The  second  approach  is  to  put  all  the  potential  h-arcs  in  first 
and  then  take  out  the  potential  h-arcs  one-by-one,  where  each  deletion 
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will  decrease  the  cost:  until  the  correct  number  of  h-arcs  is  reached. 
Any  further  deletions  will  increase  Ihe  cost.  Unfortunately,  even  if 
all  h-arcs  are  positive  with  respect  to  their  subpolygon,  the  parti¬ 
tion  may  not  be  optimum.  In  Fig.  5b,  each  h-arc  is  positive 
with  respect  to  its  local  subpolygon  but  the  partition  is  not  optimum. 
(Note  that  positiveness  of  an  h-arc  in  a  quadrilateral  is  the  same  as 
stability.  But  the  idea  of  stability  applied  to  vertical  arcs  as  well.  ) 
This  means  that  we  cannot  guarantee  an  optimum  partition  simply 
because  no  h-arc  can  be  deleted  one  at  a  time. 


fig.  5.  Counter  examples  for  the  intuitive  approach 


Let  us  outline  the  idea  of  an  O(n)  algorithm  for  finding  the  S -optimum 

partition  of  a  monotone  basic  polygon.  First,  wc  get  all  the  potential 

h-arcs  by  the  one -sweep  algorithm.  Then,  we  start  from  the  highest 

potential  h-arc  and  process  each  potential  h-arc  from  the  highest  to  the 

lowest.  For  each  potential  h-arc,  we  try  to  get  the  l  -optimum  partition 

of  the  upper  subpolygon  of  that  arc  (i.  e.  the  i  -optimum  partition  of  the 

subpolygon  bounded  below  by  that  h-arc).  The  t  -optimum  partition  in 

the  subpolygon  is  obtained  by  comparing  the  cost  of  the  f -optimal  partition 

among  the  partitions  of  the  upper  subpolygon  which  contain  one  or  more 

potential  h-arcs  with  that  of  the  fan  in  the  upper  subpolygon. 

If  we  try  all  possible  combinations  of  the  potential  h-arcs  as 

3 

candidates  for  the  1  -optimal  partitions,  we  need  0(n  )  operations  to 
find  the  i  -optimum  partition.  Fortunately,  there  are  some  dependence- 
relationships  among  these  potential  h-arcs.  Hence,  certain  subsets  of 
the  potential  h-arcs  will  either  all  exist  or  all  disappear  in  the  f  - 
optimum  partition  of  the  monotone  polygon.  We  shall  be  dealing  with 
potential  h-arcs  most  of  the  time,  so  we  shall  use  "arcs"  instead  of 
potential  h-arcs  for  brevity. 

Consider  the  monotone  basic  polygon  shown  symbolically  in 

Fig.  6.  There  arc  three  potential  h-arcs,  denoted  by  h,  ,  h.,  and  h.  . 

k  j  i 

V  is  the  local  maximum  vertex  and  V,  is  the  local  minimum  vertex, 
n  1 

Without  loss  of  generality,  we  can  assume  w  <  w  '  for  a  =  i,  j  and  k. 

a  a 

Since  we  shall  deal  with  subpolygons  bounded  by  two  potential  h-arcs, 

let  us  use  h  for  V  and  h  for  V,  (i.  e.  we  consider  these  vertices  as 

n  n  I  1 
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degenerated  arcs).  From  Lemmas  1  and  3,  the  I -optimum  partitions  of 
the  subpolygons  bounded  by  two  potential  h-arcs  (i.c.  the  white  area  of  the 
polygon  in  Fig.  6)  are  all  fans. 

Assume  (i)  h^  is  positive  in  the  subpolygon  bounded  by  h  and  h. 

but  h,  is  negative  in  the  subpolygon  bounded  by  h  and  h.  , 
k  n  i 

(ii)  h.  is  positive  in  the  subpolygon  bounded  by  h,  and  h.  but  h. 

J  k  a  j 

is  negative  in  the  subpolygon  bounded  by  h^_  and  h^  ,  and 

(iii)  h.  is  positive  in  the  subpolygon  bounded  by  h.  and  h^  only. 
Then  either  the  three  arcs  h^,  h.,  h^  all  exist  or  no  h-arcs  exists  in  the 
optimum  partition. 

This  shows  that  the  existence  of  an  h-arc  depends  on  the  existence 
of  another  h-arc. 

We  shall  use  the  notations 

to  denote  the  cost  of  the  i  -optimum  partition  of  the 
subpolygon  bounded  above  by  h.  and  bounded  below 
by  h. ,  and 

to  denote  the  cost  of  the  fan  in  the  subpolygon  bounded 
above  by  h.  and  bounded  below  by  h.  . 


Fig.  6.  An  octagon  with  three  potential  h-arcs. 


In  Fig.  6,  the  condition  for  h^  to  be  positive  with  respect  to  the 
whole  polygon  is  (compare  (5a)) 


i:) 


i  i  \  r  <  w, 

|wk:wk>  '  "V  "k  1 


(7) 


The  L.HS  of  (7)  is  denoted  by 

fh 


i:) 


and  is  called  the  supporting  weight  of  h,  with  h  as  the  ceiling  (the 

k  n 

definition  of  ceiling  will  be  given  formally  later).  Note  that  the  L.11S  of 
(7)  depends  only  on  the  weights  of  vertices  in  the  upper  subpolygon  of  h^  , 
In  terms  of  the  supporting  weights,  we  can  write  the  three 
conditions  (i),  (ii)  and  (iii)  as  follows: 
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■mtilfcwMfii 


.  ^ 


(i) 


(ii) 


(iii) 


An  arc  h.  is  a  son  of  tiie  arc  h.  (or  h.  is  the  father  of  h.) 

J  -  11  -  J 

if  the  following  conditions  are  satisfied: 

(i)  h.  is  above  h.  (the  son  is  above  the  father) 

(ii)  In  any  subpolygon  containing  h.  and  h.  ,  the  arc  h.  will  exist 
in  the  i  -optimum  partition  of  the  subpolygon  if  and  only  if  In 
exists  in  the  l  -optimum  partition. 

(iii)  In  is  the  highest  arc  that  satisfies  (i)  and  (ii). 

It  is  easy  to  see  that  every  arc  can  have  at  most  one  father  but  an  arc 
can  have  many  sons.  Also  the  ancestor-descendant  relationship  is  a 
transitive  relationship.  If  an  arc  exists  in  the  l  -optimum  partition,  all 
its  descendants  will  also  exist. 

An  arc  h^  is  a  ceiling  of  an  arc  h.  if  the  following  conditions 
are  satisfied; 

(i)  h,  is  above  h. 

k  1 

(ii)  h,  is  not  a  descendant  of  h. 

k  -  i 

(iii)  h  is  the  lowest  arc  which  satisfies  (i)  and  (ii). 

K 
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Consider  two  partitions  of  a  subpolygon  as  shown  in  Fig.  7. 


(The  shaded  areas  are  optimally  partitioned 
and  the  blank  areas  are  partitioned  by  a  fan. 
The  h-arcs  in  the  shaded  area  are  all 


descendants  of  h.  .  ) 
J 


The  cost  of  partition  of  Fig.  7a  is 


where  the  cost  of  partition  in  Fig.  7b  is 


The  condition  for  the  partition  in  Fig.  7a  to  be  cheaper  than  lhat  in  Fig.  7b 
is  (similar  to  (5a)) 


S 


<  w. 

1 
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In  order  to  give  an  intuitive  meaning  of  the  supporting  weight  S 


/\\ 
\h.  b 

'  i' 


let  us  regard  h^  and  h.  in  Fig.  7  as  fixed  while  the  position  of  h.  can  be 

moved  up  or  down  by  increasing  or  decreasing  the  values  of  w.  and  w. '. 

If  h  moves  up  and  coincides  with  h.  ,  i.e.,  w.  =  w.  ,  the  partition  in  Fig.  7a 
i  J  i  J 

costs  less  than  or  equal  to  the  partition  in  Fig.  7b.  If  the  position  of  In 

moves  down  gradually  from  h  ,  there  will  be  a  position  for  which  the  cost 

of  the  partition  in  Fig.  7a  is  equal  to  the  cost  of  the  partition  in  Fig.  7b. 

We  can  consider  this  position  as  a  fictitious  arc  f.  ,  i.e. 

- - 


cCk)tHo(fj)  ="o(fk)  ■ 


the  t  -optimum  partition  of  the  subpolygon  bounded  by  h  and  h.  becomes 

K  1 

a  fan.  The  arc  f.  is  called  the  floor  of  h.  .  Note  that  the  minimum  of  the 
J  -  J 

two  weights  associated  with  f.  is  the  supporting  weight  of  h.  . 

j  >1 

We  now  give  two  examples  to  illustrate  the  concepts,  notations 
and  the  algorithms.  Then  a  formal  description  of  the  algorithm  will  be 


given. 


Consider  a  monotone  basic  polygon  with  five  potential  h-arcs, 


h  ,h  . h  where  h,  is  the  highest  arc  as  shown  symbolically  in 

b  5  2  u 

Fig.  8  Let  w  <  w'  for  i  =  a,b . e.  The  maximum  vertex, 

p  i  l 

which  lies  above  h^  ,  has  the  weight  w^.  and  the  minimum  vertex, 
which  lies  below  h^  ,  has  the  weig  t  w  .  Wc  can  regard  w^.  (and  w^) 
as  a  degenerated  arc  and  vise  h^  to  represent  w^.  (and  h^  to  repre¬ 
sent  w  ) . 
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Fig.  8.  A  12-gon  with  5  h-arcs. 


Example  1 


Let  us  write  down  the  comparisons  made  in  the  algorithm 
First,  we  compare 


ZZ 


Fig.  9.  Illustrations  for  Kxample  1. 

9a.  To  find  f,  . 

o 

In  the  equation,  is  the  only  unknown.  In  compulation,  we  do 

not  use  the  equation  but  use  the  supporting  weight  of  instead  (h  is 

the  ceiling  of  h,  ).  If  the  h-arc  ii  .  is  below  or  coincides  with  f.  ,  which 
6  6  b  6 

means  that  is  negative  with  respect  to  the  smallest  subpolygon,  h^ 
should  be  deleted  and  never  appear  in  the  f  -optimum  partition.  For 
simplicity,  we  shall  assume  all  arcs  and  floors  have  distinct  positions 
in  the  example. 

Let  us  assume  that  f^  is  below  hr,  or  symbolically  we  write 

Vf6  ' 


Fig.  9b.  The  position  of  f  . 


Then  we  do  the  next  comparison, 


Fig.  9c.  To  find  f  . 


Assume  that  f^/fg,  i.  e.  is  a  son  of  and  h4/f5  ,  the 
comparison  is 


Fig.  9d.  Condense  h.  to  li  and  find  f 

o  S  f  \  S 


Note  that  f  is  in  a  sense  the  combined  floor  of  h  and  h  anc 
ha  1 1  '5 

becomes  the  ceiling  of  h  .  The  equation  can  also  be  written  as 

5 


If  h 4 / f 6  _  ,  the  next  comparison  will  be 


Fig.  9c.  To  fi  nd  f  ,  . 

4 


Assume  that  f,  Jl.  ,  i.e.  h  is  a  son  of  h  and  li  /f  ,  we  have 


Assume  that  f/c  i.  e.  h0  is  a  son  of  h,,and  h, /f,  ,  we  have 

oo4jo  j  £  ’  1  c. 


,h 

t  /  h  _  v 

(J 

\  +  H  /  2 

=  H  J,  7 

\h2 

'  °  '  65432 ' 

0 'f6543a' 

and  is  the  ceiling  of  h^  .  Now  if 


the  partition  consisting  of 


is  the  9  -optimum  partition. 
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Fig.  9k.  The  l  -optimum  partition. 


will  he  the  l  -optimum  partition. 


Example  Z.  The  successive  comparisons  are 


Fig.  10.  Illustrations  for  Example  Z. 
10a.  To  find  f,  . 


and  h  becomes  the  ceiling  of 


Assume  the 

1  6  5  2 


the  partition  j  is  the 


C  -optimum  partition. 


Fig.  lOh.  The  f  -optimum  partition. 

Had  we  assumed  {  /i  and  f  0/h,  then  both  h  and  h, 
52  6  52  1  5  2 

should  also  be  removed  and  we  are  left  with 


f.  against  h  . 

6  1 


If  h  /f  , 
1  6 


then  we  have  the  (  -optimum  partition 


ing  of 


From  the  above  two  examples,  we  can  see  that  h^ 

h.  if  h,  is  tli c  lowest  arc.  above  li  such  that  the 
i  k  i 


is  the  ceil 
supporting 


weight  of  lv  is  smaller  than  or  equal  to  that  of  h.  . 


Let  us  outline  the  algorithm  for  finding  the  (  -optimum  partition 
of  a  monotone  basic  polygon. 

1.  Get  all  the  potential  h-arcs  of  the  polygon  by  the  one-sweep 

algorithm.  (All  the  h-arcs  form  a  list  with  the  arc  V  -V  ,  at 

b  b 

the  bottom.  ) 

2.  Process  the  potential  h-arcs  one  by  one,  from  the  top  to  the  bottom. 
(We  try  to  find  the  l  -optimum  partition  of  the  subpolygon  bounded 
below  by  the  arc  being  processed.  ) 

2a.  Let  h  be  the  arc  currently  being  examined,  h  be  the  arc 
immediately  above  h^,  and  h^  be  the  arc  immediately  below 
h_  in  the  list.  If  hn  is  negative  with  respect  to  the  subpoly- 

K  lx 

gon  bounded  above  by  h^.  and  below  by  h^,  delete  h^,  other¬ 
wise  go  to  Step  2c. 

2b.  Once  h  and  its  descendants  are  deleted,  we  backtrack  to  h_ 

XX 

and  compare  the  cost  of  the  partition  with  h^,  and  its  descend¬ 
ants  against  the  cost  of  the  fan  in  the  subpolygon  bounded  above 
by  the  ceiling  of  h^,  and  below  by  h^.  If  the  fan  is  t  -  optimum 
in  the  subpolygon,  we  will  delete  h^.  and  repeat  this  step  until 
no  further  deletion  is  possible.  Then  we  move  to  examine  h^ . 
(The  actual  comparisons  are  done  in  terms  of  the  supporting 
weights.  ) 

2c.  Here,  h^  is  positive  in  the  smallest  subpolygon  bounded  by 
potential  h-arcs.  We  will  backtrack  to  condense  all  its 
descendants  to  h,,  as  follows.  Let  hi  be  the  ceiling  of  h„.  If 
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Wo  will 


<  SI 


'hc 


R/ 


,  h  becomes  a  son  of  h  . 

lV 


combine  h  as  well  as  all  its  descendants  to  h  and 

XV 


hc 


recalculate  the  combined  supporting  weight  S( 

\  R/ 

Replace  h^  by  h^  and  compare  the  cost  of  the  partition 
with  h  and  its  descendants  against  that  of  the  fan  in  the  sub- 
polygon  bounded  above  by  the  new  h  ,  i.  e.  h'  ,  and  below 
by  h^  .  If  the  fan  is  i  -optimum  in  the  subpolygon,  we 
delete  h^  as  well  as  its  descendants,  and  go  to  Step  2b  to 
see  if  we  can  delete  more  arcs.  Otherwise,  we  repeat  this 


step  to  see  if  we  can  condense  more  arcs. 

/h. 


2d.  Now  we  have  S| 


;:)•  <1 


,  the  supporting  weight  of 
/h. 


h^  .  The  arc  h^  is  the  ceiling  of  h^  and  S 


C 

*R , 


is  the 


supporting  weight  of  h^  .  We  move  and  process  h^. 


Before  a  formal  description  of  the  algorithm  is  given,  a  procedure 
to  process  the  list  of  potential  h-arcs  in  a  monotone  polygon  is  presented. 


Pj-f  'cdur 
In  put: 


Output: 

§te_P_S. 


Stc  p  I 


■e  MONO- PARTITION  (L) 

consists  of  a  list  of  potential  h-arcs,  passed  to  the  procedure 


via  the  argument  L.  Let  be  the  lowest  arc  in  L,  the 


one  immediately  above  h^  be  h^  .  and  h  ^  be  the  highest 


arc  in  L.  (Note  that  h.  and  h  .  are  degenerated  arcs  wi 

1  pH 

the  minimum  vertex  and  the  maximum  vertex  of  the  polygon.  ) 
consists  of  all  the  potential  h-arcs  that  exist  in  the  f -optimum 
partition  of  the  polygon. 


th 


hc  :  ~  hP+i  : 


h  r  :  =  h  ; 
K  p 


hN:=  hp-l  : 


MIN-WEIGHT  :  =  minimum  of  the  two  vertices  of  h  : 

N 


Commc nt.:  is  the  arc  to  be  processed  and  hc  is  the  ceiling 

of  the  subpolygon,  h  is  the  arc  immediately  below  h  in  L. 

IN  R 


Calculate  S 


—  S\h  /  a  min-weight 

11 7 


then  go  to  Step  Z 
else  go  to  Step  3  . 


4 
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While  (h„  i  h  ,  ,  )  And 
R  pH  - 


(the  supporting  weight  of  hR  s  MIN-WEIGHT)  Do 
Begin 

Remove  h  and  all  its  descendants  from  L; 


h„  :  =  the  ceiling  of  the  new  h 
End; 

Go  to  Step  4. 


If  (h  i  h  )  and  (the  supporting  weight  of  h  <  the 
—  C  p  f  1  -  K 


supporting  weiglit  of  h  ) 

V/ 


then 


Begin 

Condense  h^  and  all  its  descendants  into  h 
h^.  :  =  the  ceiling  of  h^.  '< 
go  to  Ste  p  1 ; 


II  ’ 


End 


else 


Begin 


Record  S 


,  .  as  the  supporting  weight  of  h  and  h  as 

hR/  lv  C 


the  ceiling  of  h  ; 


go  to  Step  4; 


Step  4 


"  'VS 


then 

Begin 


h  :  =  the  arc  immediately  below  the  new  h  ; 

N  H 

MIN-WElGHT  :  =  minimum  of  the  two  vertices  of 
the  new  h^  ; 
go  to  Step  1; 

End 

else  go  to  Step  5  ; 

Step  5  Exit  procedure  and  return  L  to  caller. 

Now  we  can  give  the  algorithm  for  finding  the  l  -optimum  parti¬ 
tion  of  a  monotone  basic  polygon. 

Algorithm  I 

Input  consists  of  n  positive  integers,  which  are  the  weights  of  the 

n  vertices  of  the  monotone  n-gon.  W[l]  is  the  weight  of  the 
minimum  vertex  and  W[i+l]  is  the  neighbor  of  W[i]  of  the 
n-gon  going  in  the  clockwise  direction.  Let  the  weight  of  the 
maximum  vertex  he  W(t], 


Output  consists  of  a  list  of  potential  h-arcs  which  will  exist  in  the 
l -optimum  partition  of  the  n-gon,  the  partitions  in  the  sub¬ 
polygons  bounded  by  every  two  consecutive  arcs  in  the  list 
are  fans. 
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Stop  0  For  i  :  =  Z  step  1  until  N  do 

i  -  1 

CP[  i]  :  =  V  W  lj]  •  W  [  j  + 1  ] ; 

j  =  1 

CP[  1]  :  =  0; 

Comment:  The  sum  of  adjacent  products  W[il:  W  [j]  can  be 
obtained  from  CP[j]  -CP[i]  for  1  <  i  <  j  “S  N  and  hence 
we  can  calculate  the  supporting  weights  easily. 

Step  1  Apply  the  one-sweep  algorithm  to  obtain  a  list  of  arcs. 

Fet  this  list  lie  L. 

Comment:  L  contains  (n-3)  arcs  which  includes  all  potential 
h-arcs  in  the  monotone  n-gon,  and  these  arcs  are  layered, 
one  above  another. 

Step  Z  From  L,  remove  those  arcs  which  are  not  potential  h-arcs; 

If  L  is  empty 
then  go  to  Step  6 
else  go  to  Step  3. 

Step  3  Let  the  lowest  arc  in  L  be  h^  ,  the  one  immediately  above 

h^  be  I-I3  ,  and  so  on: 

Let  the  highest  arc  in  L  be  h  ; 

Insert  h  with  weight  W[l]  below  h  ■ 
l  ^ 

Insert  h  with  weight  W|tl  above  h 
pt  1  p 

Comment:  h  is  the  ceiling  of  h 

-  Ptl  ^  p 
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Step  4  MONO- PARTITION  (L); 

Comment:  when  returned  from  MONO- PAR  TIT ! ON ,  L  will 
contain  all  the  ceiling  arcs  with  their  descendants  in  the  i-opt 
mum  partition. 

Step  5  Remove  h,  and  h  ,  from  L; 

- - —  1  pH 

Step  6  Output  L  and  stop. 

This  algorithm  has  been  implemented  in  Pascal  and  the  listing  of 
the  computer  program  is  given  in  Appendix  I. 

Lemma  5.  Any  arc  which  is  deleted  from  the  arc-list  L  in  Step  2  of 
the  procedure  MONO -PARTITION  cannot  be  present  in  the  i  -optimum 
partition  of  the  polygon. 

Proof.  There  are  two  cases  in  which  an  arc  is  deleted  from  L: 

(1)  Its  ancestors  are  deleted.  It  follows  from  the  definition  of  the 
ancestor -descendant  relationship  that  it  cannot  be  present  in  the 

I  -optimum  partition  of  the  polygon. 

(2)  It  is  the  h^  which  satisfies  the  logical  condition  of  Step  2  of  the 
procedure.  Hence,  in  the  subpolygon  bounded  below  by  h^  and  above  by 
h^,  the  partition  with  h^  and  its  descendants  costs  more  than  or  equal  to 
that  of  the  fan.  Hence,  the  partition  with  h  and  its  descendants  is  not 

!  -optimum  in  the  subpolygon  and  h  as  well  as  its  descendants  should  not 
appear  in  the  1  -optimum  partition  of  the  whole  polygon.  ■ 


Lemma  6.  After  an  arc  fu  has  been  processed,  the  subpolygon  be¬ 
tween  h.  and  its  ceiling  is  optimally  partitioned. 

Proof.  The  h-arcs  remaining  in  the  partition  of  the  subpolygon  are  all 
descendants  of  h.  .  By  definition  of  the  ancestor-descendant  relationship, 
the  partition  of  the  subpolygon  is  optimum.  ® 

Lemma  7.  Let  V  be  the  maximum  vertex,  and  h,  ,  h,  ,  .  .  .  ,  h.  be 
-  t  k  k- 1  jfl 

a  set  of  h-arcs  in  the  partition  such  that 


h  /h  /  . . .  /h.  / h.  , 

k  k-1  jH  J 


h  is  the  ceiling  of  h^  ^ 


h.  , .  is  the  ceiling  of  h.  , 
J  +  l  J 


then  the  supporting  weights  of  these  h-arcs  satisfy 


Proof.  Assume  that  one  of  the  i ncqualitic s  is  not  satisfied,  say 


*  -ft- 


Then  if  h.  exists  h.  will  also  exist,  h.  .  becomes  a  son  of  h  .  This 
J  J  1  *  .)  1  1  j 

contradicts  the  assumption  that  h.  .  .  is  a  ceiling  of  h.  .  ■ 

J  +  !  .1 


Lemma  8.  Any  arc  which  remains  in  L  at  the  end  of  the  procedure  must 
be  present  in  the  t  -optimum  partition  of  the  polygon. 

Proof.  We  can  divide  the  h-arcs  in  L  at  the  end  of  the  procedure  into  two 
groups: 

(i)  those  which  are  descendant  of  some  other  arcs  in  the  output,  and 

(ii)  those  which  have  no  ancestor  in  the  output. 

By  the  definition  of  the  descendant-ancestor  relationship,  the  arcs  in 
group  (i)  must  be  present  in  the  f  -optimum  partition  whenever  their  cor¬ 
responding  ancestors  in  group  (ii)  is  present  in  the  f  -optimum  partition. 
Hence,  we  have  only  to  show  that  all  arcs  in  group  (ii)  must  be  present 
in  the  l -optimum  partition. 

Let  be  the  maximum  vertex  and  the  set  of  arcs  in  group  (ii) 

be  h.  ,  h  ,,...,  h. h.  such  that  h,  /h,  ,  /  •  •  •  /h.  . ,  /h.  .  Since  none 

k  k-1  j+1  j  k  k-1  j+1  j 

of  these  arcs  has  an  ancestor,  we  must  have 

h  as  the  ceiling  of  h  , 
k  k*  1 


and  h. . ,  as  the  ceiling  of  h.  . 

J  +  1  J 

From  the  logical  condition  in  Step  1  of  the  procedure, 


we  have 

(10) 


From  Lemma  7  and  (10),  we  have 


which  implies  that 


>  C 


’*)♦<:(>  )  (Vl) 
v  Kv  °'k,  1 


Kim:;., 


+  •  •  •  f  c 


c)  -.ft)] 


=  the  cost  of  the  X  -  optimum  partition  of  the  polygon. 


In  other  words,  for  any  arc  h.  in  group  (ii)  of  L,  i  =  k,  k- 1 . j  + 1 ,  j  , 

all  the  arcs  above  h  in  L  must  be  present  in  X  -optimum  partition  of  the 

upper  subpolygon  of  h.  .  Since  (h*)*  c(h*)  *  Ho(hj)  ■thcya11 

should  be  present  in  the  X  -optimum  partition  of  the  monotone  basic 


polygon. 


Theorem  Z.  The  partition  obtained  by  the  algorithm  is  S  -optimum. 

Proof.  From  Theorems  3  and  4  of  Part  I,  we  know  that  all  the  h-arcs 
present  in  the  X  -optimum  partition  are  potential  h-arcs  and  hence  are 
included  in  the  arc -list  L  obtained  by  the  one -sweep  algorithm-  It  follows 
from  Lemmas  5  and  8  that  all  the  arcs  which  are  deleted  from  L  cannot  be 
present  in  the  X  -optimum  partition  and  all  the  arcs  which  remain  in  L  must 
be  present  in  the  X -optimum  partition.  Further,  from  Lemma  1,  the 
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I  -optimum  partition  in  any  subpolygon  bounded  by  two  adjacent  potential 
h-arcs  in  L  must  be  a  fan.  Hence,  the  partition  consisting  of  the  h-arcs 


•• 


output  by  the  algorithm  and  with  fans  in  every  subpolygons  bounded  by 
two  adjacent  arcs  in  L  must  be  i  -optimum.  ■ 

Let  us  examine  how  much  time  we  spend  in  executing  the 
algorithm. 

Step  0  and  Step  1  each  scans  the  polygon  once,  and  lienee  takes 
O(n)  time.  Since  there  are  at  most  n- 3  arcs  in  L,  Step  2  also  takes 
O(n)  time.  There  are  three  nested  loops  in  the  procedure.  The  inner¬ 
most  one  is  in  Step  6,  the  middle  one  spans  from  Step  1  to  Step  3,  while 
the  outermost  one  spans  from  Step  1  to  Step  5.  Whenever  the  innermost 
loop  is  executed  once,  a  potential  h-arc  is  deleted  from  L.  Whenever 
the  middle  loop  is  executed  once  (i,e.  the  "then"  part  of  Step  3  is  exe¬ 
cuted  once),  a  potential  h-arc  is  condensed  into  its  father.  Once  an  arc 
is  deleted  or  condensed,  it  will  never  be  examined  again.  Since  there 
are  at  most  n-3  potential  h-arcs  in  L,  the  total  number  of  executions  in 
Step  2  and  Step  3  is  O(n).  The  outermost  loop  will  also  be  executed  at 
most  (n-3)  times.  Hence  the  whole  algorithm  will  finish  its  work  in  O(n) 
time . 


3.  The  Convex  Polygon 

There  may  be  several  local  maximum  vertices  in  a  general  con¬ 
vex  polygon.  Let  us  still  draw  the  polygon  in  such  a  way  that  the  global 
minimum  vertex  is  at  the  bottom.  From  Theorem  4  of  part  I,  we  know  that 
all  potential  h-arcs  are  still  compatible  in  a  general  convex  polygon.  How¬ 
ever,  unlike  those  in  a  monotone  polygon,  the  potential  h-arcs  no  longer 
form  a  linear  list.  Instead,  they  form  a  tree,  called  an  arc-tree.  In  Fig.  11a, 

there  is  a  12-gon  with  6  potential  h-arcs  and  they  are  labelled  as  h  ,  h  h  h 

c  3  4  5 

h.  ,  and  h  .  (Note  that  we  also  obtain  V  -V  V  -V  and  V.-V  from 
b  7  4  3  7  boo 

the  one- sweep  algorithm.  In  order  to  have  a  simpler  example,  let  us 

assume  that  all  these  three  arcs  are  unstable  and  hence  are  not  shown 

in  Fig.  11a,  )  To  get  a  better  feeling  of  the  arc-tree,  we  can  redraw  the 

12-gon  as  shown  in  Fig.  lib.  Again,  we  regard  as  a  degenerated 

arc  h  ,  V  as  a  degenerated  arc  h  and  V  as  a  degenerated 
lll.  o  11 

arc  h^  . 

The  father-son  relationship  still  holds  for  the  h-arcs  in  a  gen¬ 
eral  polygon,  and  we  can  also  define  supporting  weights  of  the  arcs  in  a 
similar  way.  The  only  difference  is  that  the  ceiling  of  a  subpolygon  may 
consist  of  more  than  one  arc.  Before  we  can  calculate  the  supporting 
weight  of  any  arc,  we  must  process  all  the  arcs  above  it,  j.c.  all  the  arcs 
in  its  upper  subpolygon.  Hence,  we  can  do  a  post-order  traversal  through 
the  arc  tree.  Let  us  consider  the  following  two  examples,  \gain,  for 
simplicity,  we  assume  that  all  arcs  have  distinct  positions  in  the  examples. 
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(.b) 


A  gene ral  1  2-gon 


Example  3. 


We  first  compare 


Assume  h./f_ 
4  5 


we  compare 


Assume  h_/f.  and  wc  condense  h  into  h 

3  4  5  4  5 


Before  we  can  process  . 


we  have  to  process  h 


the  next  comparison  is; 


Assume  h.  /f_  ,  we  compare 

D  f 


We  have  h,/f,  and  f  / f .  ,  we  condense  h_  int 


Fig.  1 2f.  To  find  f. 


76  ' 


Assume  h^/f^  and  next  we  process  the  arc  ,  using  both  h^  and  h 
as  the  ceilings  of  h^  , 


H, 


V  h6 


+  H 


OV  f 


=  H 


h4'h6 


Fig.  1  2g.  To  fi nd  f ,  . 


If  12.'  the  I -optimum  partition  of  (lie  whole  polygon  consists 

of  all  six  h-arcs  h_,  h-,  h  ,  h  ,  h,,  and  h_  .  If  f  .  ,_/h,  ,  all  six 
i.  i  4  b  b  i  5  47 (>  3  Z  1 

h-nrcs  will  be  removed  and  the  /  -optimum  partition  is  a  fan. 


Example  4 


We  first  compare 


Fig.  13.  Illustrations  for  Example  4. 

13a.  To  find  fc  . 

5 


Assume  h./f, 
4  5 


and  we  compare 


Fig.  13b.  To  find  f  . 


Suppose;  »  we  remove  h..  as  well  as  its  descendants  h.  and  h  . 

rr  43?,  1  L  3  '1 

Assume  f.  /fr  and  f.  /h  ,  we  remove  h.  from  the  polygon.  Now,  we 
6  5  6  1  ''  ' 

have  f  /f  and  f  /h  ,  so  we  remove  h  from  the  polygon.  Finally, 

7  5  7  1  * 

we  have  h  /f  ,  and  the  f  -optimum  partition  of  the  polygon  consists  of  one 

1  5 


h-arc  hj.  . 


From  the  above  two  examples,  we  have  the  following  observa¬ 
tions. 

(1)  Before  we  can  process  a  potential  h-arc,  we  have  to  process  all 
the  arcs  above  it.  Hence,  wc  should  do  a  post-order  traversal,  starting 
at  the  root  of  the  arc  tree,  i.  e.  the  degenerated  arc  h^  . 

(2)  Whenever  wc  do  a  condensation  or  deletion,  we  always  pick  the 
ceiling  arc  which  has  the  highest  floor  first,  i.  e.  the  one  with  the 
largest  supporting  weight.  Hence,  wc  should  keep  track  of  the  order 


of  the  ceiling  arcs. 


(3)  Once  a  ceiling  arc  h.  of  h.  is  removed  or  condensed,  the  ceiling 
arcs  of  h.  become  the  ceiling  arcs  of  h.  and  we  have  to  update  the 
order  of  all  the  ceiling  arcs  of  h,  , 

One  way  of  keeping  track  of  the  order  of  the  ceiling  arcs  is  to 
keep  them  in  a  priority  queue. 

Now,  let  us  outline  the  algorithm  for  finding  the  optimum  parti¬ 
tion  of  a  general  convex  polygon. 

1.  Get  all  the  potential  h-arcs  of  the  polygon  by  the  one-sweep 
algorithm.  (All  the  h-arcs  form  a  tree,  with  the  root  at  the 
bottom.  Let  the  arc-tree  be  T.  ) 

2.  Process  the  h-arcs,  one  by  one,  from  the  leaves  to  the  root.  (We 
always  process  the  children  before  we  process  the  father,  and  we 
always  obtain  the  optimum  partition  of  the  subpolygon  bounded  below 
by  the  arc  being  processed. ) 


Let  h  be  the  arc  currently  being  examined,  Uu  be  the  set  of  arcs 
K  w 


R 


immediately  above  h  ,  and  h  be  the  arc  immediately  below  h 

Jt\  IN  K 

in  T.  If  h  is  negative  in  the  subpolygon  bounded  above  by  the 
rv 

arcs  in  U  and  below  by  h  ,  delete  h  ,  else  go  to  step  5. 

K  In  K 

4.  Once  h  and  its  descendants  are  deleted,  we  examine  the  arcs  in 

X\ 

u  to  sec  if  we  can  delete  more  arcs,  if  yes,  we  delete  tin'  arc 
with  the  largest  supporting  weighl ;  then  we  include  its  ruling  arcs 
into  U  and  repeat  this  step.  Otherwise,  we  move  to  process  the 
next  arc. 
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Now,  h„  is  positive  in  the  smallest  subpolygon.  If  there  exists 
R 

some  arc  in  ,  say  h. ,  such  that 

R  j 


<  the  supporting  weight  of  h. 


we  will  pick  the  arc  with  the  largest  supporting  weight  in 

U  ,  condense  it  with  its  descendants  into  h  and  include  all 

its  ceiling  arcs  into  U^.  Then  we  compare  the  cost  of  the  partition 

with  h^  and  its  descendants  against  that  of  the  fan  in  the  subpolygon 

bounded  above  by  the  arcs  in  U_,  and  below  by  h„,.  If  the  fan  is 

t  -optimum  in  the  subpolygon,  we  remove  h  as  well  as  all  its 

K 

descendants  from  T,  and  we  examine  the  arcs  in  U  to  see  if  we 

K 

can  delete  any  more  arcs.  Otherwise,  wc  examine  the  arcs  in  LT 

R 

to  see  if  we  can  condense  any  more  arcs. 


Now,  S 


U 


s  the  supporting  weight  of  every  arc 

U. 


arcs  in  U0  are  the  ceiling  arcs  of  li  and  S 
R  R 


R 

'R 


in  U  .  The 
R 

is  the  supporting 


R  ' 


weight  of  h 


We  move  to  process  the  next  arc. 


Before  presenting  the  algorithm,  let  us  describe  a  recursive 


procedure  to  process  the  potential  h-arcs  of  any  subpolygon. 

Procedure  PARTITION  (ROOT) 

Input:  consists  of  a  set  of  potential  h-arcs  of  a  subpolygon.  These 

arcs  are  arranged  in  the  form  of  an  arc  tree,  like  the  one 
shown  in  Fig,  lib.  The  root  of  the  tree  is  passed  to  the 
procedure  via  the  argument  ROOT. 

Output:  consists  of  a  set  of  the  potential  h-arcs  which  appear  in  the 

f  -optimum  partition  of  the  subpolygon.  We  can  divide  that 
arcs  into  two  types:  (i)  those  arcs  which  are  descendants  of 
some  other  arcs  in  the  set  and  (ii)  those  arcs  which  have 
no  ancestor  in  the  set.  The  arcs  in  type  (i)  are  con¬ 
densed  into  their  ancestors  and  can  be  traced  out  from  the 
arcs  in  type  (ii).  The  arcs  in  type  (ii)  are  called  ceiling 
arcs  and  are  kept  in  a  reduced  arc  tree.  The  root  of  the  arc 
tree  is  passed  back  to  the  caller  via  the  parameter  ROOT. 
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Step  0  Let  the  arc  at  the  root  of  the  input  arc  tree?  be?  h  ; 

M1N-WK1GHT  :  =  the  weight  of  the  minimum  of  the  two 
vertices  of  h^  ; 

T  :  =  an  arc  tree  with  only  one  arc,  ; 


Step  1  For  each  arc  immediately  above  h^  in  the  input  arc-tree  Do 
Begin 


Step  la 


Step  lb 


Step  lc 


Let  the  arc  to  be  processed  be  h  ; 

If  there  exists  a  non-dcgencrated  arc  above  h^ 
then  go  to  Step  lb 
else  go  to  Step  If; 

Comment:  h„  is  immediately  above  h  ,  . 

PARTITION  (h  )  ; 

Let  the  subtree  returned  be  T  '  ; 


Comment:  Before  processing  h,,  ,  the  subtrees  of  h.. 

K  R 

are  first  processed  recursively. 


Let  U  be  the  set  of  arcs  immediately  above  h  in  7-/  ; 
K  R 


Calculate  ; 


MIN-WKICI  IT 


then  go  to  Step  Id 
else  go  to  Stc p  l e . 


ft 


Step  Id 


Remove  h  from  T '  ; 

r\ 

while  (there  exists  a  non-degeneratecl  arc,  h.  ,  in  U.,)  and 
-  j  R  - - 

(the  supporting  weight  of  S  MIN-WK1GHT)  Do 
Begin 

Remove  h.  from  U_  ; 

J  R 

Remove  h.  from  T'; 

J 

Include  all  ceiling  arcs  of  h.  into  U  ; 

J  R 


Insert  the  forest  T'  into  T  such  that  all  arcs  in  U  are 

K. 


immediately  above  h„  in  T. 


Go  to  Step  li. 


Step  le  If  (there  exists  a  non-degenerated  arc  in  U^)  and  (its  sup¬ 
porting  weight  >  the  supporting  weight  of  h  ) 


Begin 


Among  all  the  arcs  in  U  ,  pick  the  one  with  maximum 

tv 

supporting  weight; 

Let  it  be  h.  ; 

J 

Condense  h.  into  h  and  remove  it  from  T '  ; 

J  K 

Include  all  ceiling  arcs  of  h.  into  U  ; 

j  R 

Fix  up  the  tree  T  '  so  that  all  the  ceiling  ares  of  lu  an 
immediately  above  h  in  T'  ; 
go  to  Step  1c; 


^  J 
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else 


Step  lg 


Regin 


Record  SI  I  as  the  supporting  weight  of  h  and  all 
\hR/  R 

arcs  in  U  as  the  ceiling  arcs  of  hD ;  insert  T  '  into  T 
R  K 

so  that  h  is  immediately  above  h  in  T; 

R  IN 


go  to  Step  li  ; 


Step  If  Let  h^.  be  the  degene  rated  arc  above  h„ 
^  R 


Calculate  S 


If  S^  j  a  MIN- WEIGHT 

then  go  to  Step  lg 
else  go  to  Step  lh  . 


Remove  h  ; 

R 

Insert  h ^  immediately  above  h  in  T. 


go  to  Step  li  . 


Step  lh  Record  Si  J  as  the  supporting  weight  of  h  and  h_,  as 


Step  li 


the  ceiling  arc  of  h  ;  insert  Hie  subtree  with  h  and  h 

R  K  O 

into  T  so  that  hR  is  immediately  above  h^  in  T. 


Step  2.  Return  T  with  root  stored  in  ROOT  to  caller. 


63 


Now,  the  details  of  the  algorithm  to  find  an  optimum  partition  of 


a  convex  polygon  is  presented. 

Algorithm  II 

Input  consists  of  n  positive  integers,  which  are  the  weights  of  the 

n  vortices  of  an  n -gon.  W[l]  is  the  weight  of  the  minimum 
vertex  and  W[i  +  l]  is  the  neighbor  of  W[i]  of  the  n-gon  going 
in  the  clockwise  direction. 

Output  consists  of  a  tree  of  potential  h-arcs  which  exist  in  the 
l  -optimum  partition  of  the  n-gaon. 


Step  0  For  i  :  =  2  step  1  until  N  do 

i  -  1 

CP[i]  :  =  X  WU3  •  W[j  +  1]  ; 

j  =  1 

CP[l]  :  -  0  ; 

Comment:  The  sum  of  adjacent  products  W[i]:W(j]  can  be 

obtained  from  CP[j]-  CP[i]  for  1  <  i  <  j  <  N. 


Step  1 


Apply  the  one-sweep  algorithm  to  obtain  a  tree  of  arcs.  Let 
this  arc  tree  be  T. 

Comine nt:  T  contains  all  potential  h-arcs  in  the  n-gon. 


Stc p  2  From  T,  remove  those  arcs  which  arc  not  potential 


St  ('ll  1 


Insert  the  flcj’i:  tic  ra  t  c<l  a  rc  h  with  weight  W|  l]  to  the 
bottom  of  tho  tree,  as  the  root  of  tl  1* •  fret'; 

Insert  a  degenerated  arc:  with  the:  local  maximum  weight 
at  the  tip  of  each  corresponding  branch  of  the  arc  tree. 

Step  4  PARTITION  (h  ); 

Comment:  is  the  root  of  T;  when  returned  from 

PARTITION,  T  will  contain  all  the  ceiling  arcs  with  their 
descendants  in  the  l  -optimum  partition. 

Ste  p  5  Remove  all  degenerated  arcs. 

Step  6  Output  T  and  stop. 

This  algorithm  has  been  implemented  in  Pascal  and  the  list¬ 
ing  of  the  computer  program  is  given  in  Appendix  II. 

Theorem  3.  The  partition  of  the  general  convex  n-gon  obtained  by  the 
algorithm  is  t  -optimum. 

Proof.  Using  arguments  similar  to  those  in  Theorem  2,  we  can  first 
prove  that  all  the  potential  h-arcs  which  are  deleted  from  the  arc-tree 
cannot  be  present  in  the  l  -optimum  partition,  then  we  prove  that  any  arc 
which  is  left  in  the  arc-tree  at  the  end  of  the  algorithm  must  be  present 
in  the  i  -optimum  partition.  Hence,  the  partition  consisting  of  the  h-arcs 
output  by  the  algorithm  and  with  fans  in  tho  subpoiygons  bounded  by  a 
potential  h-arc  and  the  arcs  immediately  above  it  in  the  output  arc- 
Ireo  must  be  l -optimum.  H 
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Let  ns  examine  how  much  time  we  spend  in  executing  the  algorithm. 

Steps  0  and  1  each  scans  the  polygon  once,  and  hence  takes  O(n) 
time.  Since  there  are  at  most  n-3  arcs  in  T,  Step  Z  also  takes  O(n)  time. 
There  will  be  a  recursive  procedure  call  for  each  arc  in  T  (except  the 
leaf  nodes).  Inside  each  procedure  call,  there  are  two  nested  loops.  The 
innermost  loop  is  the  "while"  loop  in  Step  Id  and  the  outer  one  spans  from 
Steps  lc  to  le.  Whenever  the  innermost  loop  is  executed  once,  a  potential 
h-arc  is  deleted  from  T.  Whenever  the  outer  loop  is  executed  once  (i.e. 
the  "then"  part  of  Step  le),  a  potential  h-arc  is  condensed  into  its  father. 
Once  an  arc  is  deleted  or  condensed,  it  will  never  be  examined  again.  In 
order  to  carry  out  the  deletion  and  conde  ation  efficiently,  we  cannot 
examine  all  the  arcs  in  UR  each  time  we  go  through  the  loop.  Hence,  we 
need  to  order  the  arcs  in  UR  in  a  priority  queue  and  it  takes  0(log  n)  to 
update  the  queue  each  time.  Hence,  it  takes  0(n  log  n)  time  in  executing 
Step  4  of  the  algorithm.  Steps  5  and  6  each  takes  O(n)  time.  Hence,  the 
whole  algorithm  takes  0(n  log  n)  time  to  find  the  t  -optimum  partition. 
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4.  A  closer  look  at  the  optimum  partitions 

Wc  now  present  some  theorems  which  enable  the  algorithm  to 
divide  the  polygon  into  several  subpolygons  and  hence  can  improve  the 
average  performance  of  the  algorithm.  These  theorems  have  also  been 
mentioned  in  [4]  without  detailed  proofs. 

Let  us  consider  the  polygons  where  there  are  two  or  more 
vertices  with  equal  weights  w  . 

Lemma  9.  For  every  choice  of  V, ,  V,,  ..  .  (as  prescribed  in  Part  I),  if 
the  weights  of  the  vertices  satisfy  the  condition 

w,  =  w,  <  w„  <  •  •  •  <  w  , 

12  3  n 

then  V  -V  exists  in  every  optimum  partition  of  the  n-gon. 

1  u 

Proof.  The  lemma  is  true  if  V,-V,,  is  a  side  of  the  n-gon.  Hence,  we 
can  assume  that  V^,  are  not  adjacent  to  the  same  side  of  the  n-gon. 

The  proof  is  by  induction  on  the  size  of  the  n-gon.  The  lemma  is 
true  for  a  triangle  and  a  quadrilateral.  Assume  that  the  lemma  is  true  for 
all  k-gons  (3  £  k  £  n-1)  and  consider  the  optimum  partitions  of  an  n-gon. 

By  Lemma  3  of  Part  I,  we  know  that  there  are  at  least  two  vertices 
with  degree  two  in  each  optimum  partition  of  the  n-gon.  We  have  the 
following  two  cases. 

(i)  In  an  optimum  partition  of  an  n-gon,  one  of  the  vertices  with  degree 
two,  say  V.,  has  weights  larger  than  w^.  In  this  case,  we  can  form  an 
(n-l)-gon  by  removing  V.  with  its  two  sides.  By  induction  assumption, 

^ 2  Present  in  every  optimum  partition  of  the  (n-l)-gon. 


(ii)  Consider  the  complementary  case  of  (i),  i.e.  all  vertices  with 
degree  two  have  weights  equal  to  w^  in  an  optimum  partition  of  the  n-gon. 

In  other  words,  V  and  V-  are  the  only  two  vertices  with  degree  two  in 
that  optimum  partition,  as  shown  symbolically  in  Fig.  14a.  Note  that 
every  arc  in  the  optimum  partition  must  dissect  the  n-gon  into  two  sub¬ 
polygons  in  such  a  way  that  V^,  can  never  appear  in  any  subpolygon 
together,  else  there  will  be  more  than  two  vertices  with  degree  two  in 
the  optimum  partition.  In  Fig.  14b  we  show  a  partition  of  the  n-gon  in 
which  V  and  V  are  connected.  Let  us  denote  the  n- 2  triangles  in  Fig.  I4a 

1  Ct 

by  P, ,  P-,  .  .  .  ,  P  Except  P,  and  P  _,  all  the  other  n-4  triangles  are 

made  up  of  one  side  and  two  arcs  each.  For  each  of  these  n-4  triangles, 
we  can  find  a  unique  triangle  in  Fig.  14b  such  that  they  both  consist  of 
the  same  side.  We  use  P /  to  denote  the  image  of  P.  in  Fig.  14b.  The 
only  two  triangles  left  unmatched  in  Fig.  14b  are  V  V  V  and  V  V_V. 
and  they  are  the  images  of  P^  and  P^  respectively.  Let  the  cost  of  P. 
be  C.  and  the  cost  of  P.r  be  C/  .  Since  C/  <  C.  for  1  s  i  <  n-2,  the 

l  ii  ii 

partition  in  Fig.  14b  is  cheaper  than  that  in  Fig.  14a  and  we  have 

contradiction.  ■ 
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Theorem  4.  For  every  choice  of  V  ,  V 2,  ...  (as  prescribed  in  Part  I),  if 
the  weights  of  the  vertices  satisfy  the  condition 

1  2  3  4  n 


then  every  optimum  partition  of  the  n-gon  must  contain  a  triangle  V^V^V 
for  some  vertex  with  weight  equal  to  w^.  Note  that  if  w^  =  w^  <  w^ 
<w^  s  .  •  .  s  w^  ,  then  every  optimum  partition  must  contain  the  triangle 
since  there  is  a  unique  choice  of  . 


P 


Proof.  Similar  to  Lemma  9.  wc  can  prove  thin  theorem  by  induction  on  the 
size  of  the  n-gon.  The  theorem  is  true  for  any  triangle  or  quadrilateral 
satisfying  the  above  condition.  Assume  the  theorem  is  true  for  all  k-gons 
(3  £  k  <  n-1)  and  consider  the  optimum  partitions  of  an  n-gon. 


From  Lemma  9,  we  know  that  V^,  are  always  connected  in  every 
optimum  partition.  Hence,  without  loss  of  generality,  we  can  assume  V^, 
to  be  adjacent  to  the  same  side  of  the  n-gon.  Again,  we  have  the  following 
two  cases. 


(i)  In  an  optimum  partition,  one  of  the  vertices  with  degree  two,  say  V  , 


has  weight  larger  than  w^.  In  this  case,  we  can  remove  with  its  sides  and 


form  an  (n-l)-gon.  By  induction  assumption,  every  optimum  partition  of  the 
(n-l)-gon  contains  a  triangle  V  V  V  for  some  vertex  V  where  w  =w  . 

1  t*  p  p  P  J 

(ii)  Consider  the  complementary  case  of  (i),  in  an  optimum  partition  of 


the  n-gon,  all  vertices  with  degree  two  have  weights  less  than  or  equal  to  w^  . 


Since  V  -V  is  a  side  of  the  n-gon,  for  n  ^  4,  either  V  or  V  (but  not  both) 

\  U  \  L. 


can  have  degree  two.  We  have  the  following  two  subcases: 


(a)  If  there  are  more  than  one  vertex  whose  weight  equals  w^,  we 


can  form  an  (n-l)-gon  by  removing  one  of  those  degree  two  vertices  whose 


weight  equals  w^.  By  induction  assumption,  every  optimum  partition  of  the 


(n-l)-gon  contains  a  triangle  V,V_V  for  some  vertex  V  with  w  -  w,  , 

1  2  p  p  p  3 


(b)  There  exists  only  one  vertex  of  weight  w y  In  this  case,  there 


must  be  only  two  vertices  with  degree  two  in  the  optimum  partition  of  the 
n-gon.  These  two  vertices  are  and  either  or  V^.  Without  loss  of 


generality,  we  can  assume  has  degree  2.  The  situation  is  shown  symboli¬ 


cally  in  Fig.  15a.  Again,  every  arc  in  the  optimum  partition  must  dissect  the 


n-gon  in  such  a  way  that  and  can  never  appear  in  any  subpolygon  to¬ 


gether.  In  Fig.  15b,  we  show  a  partition  containing  the  triangle  V  V  V  . 

1  “  J 


Using  arguments  similar  to  those  in  the  proof  of  Lemma  9,  we  can  show  that 
the  partition  in  Fig.  15b  is  cheaper  and  we  obtain  a  contradiction.  ■ 


Fig.  15 


Theorem  5.  For  every  choice  of  V^,  V^.  .  . .  (as  prescribed  in  Part  I),  if  the 
weights  of  the  vertices  of  the  n-gon  satisfy  the  following  condition. 


Wi  =  w2 


Wk< 


for  some  k,  3  5  k  <  n,  then  every  optimum  partition  of  the  n-gon  contains 
the  k-gon  V  -V^- • •  •  -  V^. 

Proof.  The  proof  is  by  induction  on  tt.  of  the  n-gon.  The  theorem  is 

true  for  any  triangle  and  quadrilateral.  Suppose  the  theorem  is  true  for  all 
polygons  with  (n-1)  sides  or  less  and  consider  the  optimum  partitions  of  an 
n-gon. 

From  Lemma  3  of  Part  I,  there  exist  at  least  two  vertices  having 
degree  two  in  every  optimum  partition.  We  have  the  following  two  cases. 


(i)  In  an  optimum  partition  of  the  n-gon,  one  of  the  vertices  with  degree 

two,  say  V.  ,  has  weight  larger  than  w^  In  this  case,  we  can  remove  the 

vertex  V  with  its  two  sides  and  obtain  an  (n-l)-gon.  By  induction  assumption, 

1 

every  optimum  partition  of  the  (n-l)-gon  contains  the  k-gon  V^-V^-  .  .  .  -V^  • 

(ii)  Consider  the  complementary  case  of  (i),  i.e.  ,  all  the  vertices  with 

degree  two  have  weights  equal  to  w^  in  an  optimum  partition.  Let  two  of 

these  vertices  be  V.,  V..  We  have  the  following  two  subcases: 

1  3 

(a)  k>  3.  We  first  form  an  (n-l)-gon  by  removing  V.  and  its  two  sides. 
There  are  (k-1)  vertices  with  weights  equal  to  w^  in  the  (n-l)-gon.  By  induc¬ 
tion  assumption,  every  optimum  partition  of  the  (n-l)-gon  contains  the  (k-1)- 

gon  which  includes  V.  as  one  of  its  vertices.  Since  V.  has  degree  two  in  the 

J  3 

optimum  partition,  its  two  neighboring  vertices,  say  V  and  V  ,  must  also 

x  y 

have  weights  equal  to  w^  and  the  arc  V^-V^  exists  in  the  optimum  partition 

(Fig.  16).  Similarly,  we  can  remove  the  vertex  V.  with  its  two  sides  V.-V 

J  3  x 

and  V.-V  and  form  an  (n-l)-gon.  By  induction  assumption,  every  optimum 

3  y 

partition  of  the  (n-l)-gon  contains  the  (k-l)-gon  formed  by  the  (k-1)  vertices 

with  weights  equal  to  w^  in  the  (n-l)-gon  and  V.  is  one  of  the  vertices  in  the 

(k-l)-gon.  Now,  by  pasting  the  triangle  V  V.V  and  the  (k- l)-gon  together, 

x  3  Y 

we  form  a  k-gon  which  includes  all  the  vertices  with  weight  equal  to  w^  in 
the  n-gon  and  this  k-gon  is  contained  in  the  optimum  partition  of  the  n-gon. 
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Fig.  16 


(b)  k  =  3.  In  this  case,  we  have  w  =  w,  =  w,  <  w .  ^  •  *  •  £  w  . 
Without  loss  of  generality,  we  can  assume  and  V ^  both  have  degree  two  in 
an  optimum  partition.  Again,  we  can  form  an  (n-l)-gon  by  removing  and 
its  two  sides.  By  Lemma  9,  V  and  V  are  connected  in  every  optimum 
partition  of  the  (n-l)-gon.  Since  V  has  degree  two,  V  -V  must  be  a  side 

6  4  J 

of  the  n-gon.  Next,  we  can  remove  with  its  two  sides  and  form  an  (n-1)- 
gon.  By  Lemma  9,  V^,  are  connected  by  a  side  of  the  n-gon.  The  situa¬ 
tion  is  shown  in  Fig.  17a.  Then,  the  partition  in  Fig.  17b  is  cheaper  because 


123 


+  T 


12y 


T13x  + 


23y 


and 


C(w,  .  w . .  )  <  C(w,,  w . .  )  .  ■ 

»  x  y  3  x  y  * 
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Fig.  17 


Now,  whenever  we  have  three  or  more  vertices  with  weights  equal 
to  w^  in  the  n-gon,  we  can  decompose  the  n-gon  into  subpolygons  by  forming 
the  k-gon  in  Theorem  5.  The  partition  of  the  k-gon  can  be  arbitrary,  since 
all  vertices  of  the  k-gon  are  of  equal  weight.  For  any  subpolygon  with  two 
vertices  of  weights  equal  to  w^  we  can  always  apply  Theorem  4  and  decom¬ 
pose  the  subpolygon  into  smaller  subpolygons.  Hence,  we  have  only  to 
consider  the  polygons  with  a  unique  choice  of  V^,  i.e.  ,  each  polygon  has 
only  one  vertex  with  weight  equal  to  w^. 

Because  of  Theorems  4  and  5,  Theorems  1  and  3  of  Part  I  can  be 
generalized  as  follows. 

Theorem  6.  For  every  choice  of  V  ,  V  , . . .  (as  prescribed  in  Part  I),  if  the 

1  Ct 

weights  of  the  vertices  satisfy  the  condition 


w  <w_<w,s*--<w  , 

l  c  3  n 


then  V^-V^  am*  Vj-V  ex*8*  *n  every  optimum  partition  of  the  n-gon. 


Theorem  7.  Let  and  be  two  arbitrary  vertices  which  are  not  adjacent 

in  a  polygon,  and  V  be  the  smallest  vertex  from  V  to  V  in  the  clockwise 
w  x  z 

manner  (V  /  V  ,  V  ^  V  ),  and  V  be  the  smallest  vertex  from  V  to  V 
w  x  w  z  y  z  x 

in  the  clockwise  manner  (V  ^  V  ,  V  ^  V  ).  This  is  shown  in  Fig.  18 

y  x  y  z 

where  we  assume  that  V  <  V  and  V  <  V  .  The  necessary  condition 

x  z  y  w 

for  V^-V^  to  exist  as  an  h-arc  in  any  optimum  partition  is 

w  <  w  s  w  <  w  ■ 

y  x  z  w 


......  Vw 

s 

V 
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Fig.  18 


From  Theorem  7,  we  know  that  any  arc  which  exists  as  an  h-arc  in 
some  optimum  partition  must  be  a  potential  h-arc.  In  other  words,  the 
h-arcs  in  every  optimum  partition  will  be  generated  by  the  one-sweep 
algorithm.  Hence,  by  modifying  the  condition  in  steps  lc  and  Id  of  the 
procedure  Partition  to  favor  partitions  with  more  h-arcs,  we  can  obtain 
other  optimum  partitions  which  consist  of  more  h-arcs  than  the  I -optimum 


partition. 
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5.  Conclusion 

The  problem  to  find  the  optimum  order  of  computing  a  chain  of 
matrices  has  been  around  for  several  years  [2]  .  It  has  been  used  as  a 
typical  example  to  illustrate  the  dynamic  programming  technique  in  many 
textbooks  [l][3]  .  In  this  paper,  a  new  approach  is  used  to  solve  the 
problem.  Instead  of  tackling  the  matrix  chain  product  problem  directly, 
it  is  transformed  into  the  problem  of  partitioning  a  convex  polygon  and  a 
tailor-made  algorithm  for  finding  the  optimum  partition  is  developed. 

The  algorithm  takes  0(n  log  n)  time  and  O(n)  space.  For  those  who  want 
to  trade  optimum  solution  for  shorter  execution  time,  an  O(n)  heuristic 
algorithm  has  been  presented  in  [5].  This  heuristic  algorithm  is  very 
simple  to  implement  and  its  error  bound  given  explicitly  as  a  function 
of  the  number  of  sides  of  the  convex  polygon  and  the  ratio  of  the 
weights  of  the  largest  vertex  to  that  of  the  smallest  vertex.  The 
worst  error  ratio  is  less  than  15%. 
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Appendix  I 


PROGRAM  OPTIMUM  ALGORITHM  FOR  A  MONOTONE  BASIC  POLYGON; 


CONST  MAX_S IZE  =  127; 

TYPE  P0S_1NTEGER 

LIST_PTR 
LIST  ELEMENT 


0  ..  32767;  (limited 
size  of 

*  LIST  ELEMENT; 

PACKED  RECORD 

HEAD,  TAIL  : 

HEAD_SMALL  : 

SUP_WE IGHT , 

COST  , 

BASE_PRODUCT , 
SIDE_PRODUCT  : 

DESCENDANT,  NEXT  : 
END; 


by  the  word- 
the  computer} 


POS_I NTEGER ; 
BOOLEAN; 


INTEGER; 
LIST  PTR 


VAR 


W,  CP  ; 
LIST,  LEAF: 
N  : 


ARRAY  ( 1  .  .MAX_SI'/,Ej  OF  INTEGER; 
LISTJPTR; 

POS  INTEGER; 


SEGMENT  PROCEDURE  INITIALIZING; 

(*  Handles  the  inputs  and  initializing  all  the  global  *) 

(*  variables,  *) 

VAR  I  :  INTEGER; 


BEGIN 

WRITELN  ('a 
WRITELN  (' 
WRITELN; 


linear  algorithm  to  find  all  the  h-arcs  in', 

'the  optimum'); 

partition  of  a  monotone  basic  polygon', 

'  (7/2/80)'); 


(obtain  the  inputs} 

WRITE  ('Please  enter  the  size  of  the  polygon  (between  3', 

'  and  '  ,MAX_SJ ZE-1  ,  '  )  :  '); 

READLN  (N) ; 

WR I  TEI.N ; 

WRITELN  ('Now,  starting  from  the  smallest  vertex  and  in', 

'  the  '  )  ; 

WRITELN  ('  clockwise  direction,  enter  the  weights  of', 

'  the  ver  t ices  :  ' ) ; 

FOR  I  :=  1  TO  N  DO  READ  (W(lj); 

READLN; 

WRITELN; 
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{calculate  the  cumulative  aujacenl 

products  around  the  polygon} 

Cl*  l  1  I  :=  U; 

FOR  1  :=  2  TO  N  DO  CP [ I }  :=  CP(I-1J  +  W [ 1 - 1 }  *  W [  I  ]  ; 

{initialize  the  psuedo  arc} 

NEW  (LEAF)  ; 

WITH  LEAF"'  DO 
BEGIN 

BASE_PRODUCT  :=  0; 

SIDE_PRODUCT  :=  0; 

END; 

{set  up  the  output  headings} 

WRITELN  ('the  potential  h-arcs  in  the  partitions  are  :  '); 
END;  {initializing} 


SEGMENT  PROCEDURE  0NE_SWEEP  (VAR  L  :  LIST_PTR) ; 

(*  Sweep  the  polygon  once,  collects  all  potential  h-arcs,  *) 
(*  puts  them  in  a  list.  The  address  of  the  head  of  the  *) 
(*  list  is  stored  in  L.  *) 


VAR  STACK 


TOP_ELEMENT,  SEC0ND_EL.EMKNT , 
CURRENT,  TOS 
P,  ARC  LIST 


:  ARRAY  [ 1 . . MAX_S I ZE)  OF 

POS_INTEGER; 

;  POS  INTEGER; 

:  LIST  PTR ; 


PROCEDURE  PUSH  (C  :  INTEGER)  ; 

(*  Pushes  the  index  C  onto  the  stack  and  updates  the  *) 

(*  variables  TOS,  TOP_ELEMF.NT ,  and  SECOND_ELKMENT .  *) 

j********************************************************) 
BEGIN 

STACK [TOS J  :=  C; 

SECOND_ELEMENT  :=  TOP_ELEMLNT; 

TOP_F.  LEM  ENT  :=  C; 

TOS  TOS  -  1; 

END;  {push} 


P  ROC  E  DU  R  E  PO  P  _  S  T  AC  K  ; 

^*****Jc******JT*;*******:*******:***  A  **************  **********) 

(*  Pops  the  top  element  off  the  stack  and  update:',  the  *) 
(*  variables  TOS ,  TOP_E  LEMENT ,  and  SECOND  E  LEMENT .  *) 

J********************************************************) 

BEGIN 

TOS  :=  TOS  +  1; 

TOP_ELEMENT  :=  SECOND_ELEMENT ; 

SECOND_E  LEMENT  :=  STACK[TOS  +  2); 

END;  {pop  stack} 


(*  One-sweep  begins  here.  *) 

I**********************************************************) 
BEGIN 

{initialize  the  local  variables) 

TOP  ELEMENT  :=  0; 

SECOND_ELEMENT  :=  0; 

STACK [N+l]  :=  0; 

TOS  : =  N ; 

ARC_LIST  :=  NIL; 

PUSH  ( 1 ) ; 

PUSH  (2)  ; 

CURRENT  :=  3; 

{scan  through  the  polygon  in  the  clockwise  direction} 

WHILE  CURRENT  <  N  DO 

IF  ( W  [SECOND_ELEMENT]  <=  WlTOP_ELF.MF.NT]  )  AND 

( W  l TOP_ELEMENT]  >  W [CURRENT]  ) 

THEN 
BEGIN 
NEW ( P ) ; 

WITH  P~  DO 
BEGIN 

HEAD  :=  S  ECOND_E  LEMENT ; 

TAI  I.  :=  CURRENT; 

HEAD_SMALL  :=  W (HEAD)  <=  W(TAIL) ; 

BAS E_PR0DUCT  :=  W [ H EA D ]  *  W[TAIL]  ; 

S IDE_PRODUCT  CP(TAIL)  -  C P [H EAD }  ; 

DESCENDANT  :=  NIL; 

NEXT  :=  ARC  LIST; 

END; 

A  RC_L I ST  :=  P; 

POP  STACK; 

IF  TOS  >=  (N - 1 )  {there  are  less  than  2 

elements  on  the  stack) 


1 


THEN 

BEGIN 

PUSH  (CURRENT)  ; 

CURRENT  :=  CURRENT  +  1; 

END; 

END 

ELSE 

BEGIN 

PUSH  (CURRENT) ; 

CURRENT  :=  CURRENT  +  1; 

END; 

WHILE  (TOS  <=  (N-3 )  ) 

AND  (W  [  SECOND_ELEMENT]  <=  W  ( TOP_E  LEMEN'i’l  ) 

AND  ( W[TOP_ELEMENT)  >  W[N]  )  DO 

BEGIN 
NEW ( P)  ; 

WITH  P~  DO 
BEGIN 

HEAD  :=  S  ECOND_E  LEM ENT ; 

TAIL  :=  N; 

HEAD_SMALL  :=  W  ( F1EAD]  <=  W(TAIL); 

BAS E_PRODUCT  :=  W(HEAD)  *  W[TAIL); 

S IDE_P RODUCT  :=  CP(TAIL]  -  C P [H EAD]  ; 

DESCENDANT  :  =  NIL; 

NEXT  :=  A  RC_L I ST ; 

END; 

A  RC_LIST  :=  P; 

POP_STACK ; 

END; 


L  :=  A  RC_LIST  ; 
END;  {one_swccp} 
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(VAR  L  :  LIST  PTK)  ; 

**********“3F************' 


SEGMENT  PROCEDURE  MONO^PARTITION 

(*  Obtains  all  the-  h-arcs  in  the  <’>|  >  t  i  in  t  im  partition  <>f  tin 
(*  polygon  and  returns  them  in  a  list.  The  address  of 
(*  the  head  of  the  list  is  stored  in  L.  *) 


*) 

*) 


FUNCTION  FAN_COST  (HR,  HC  :  LIST_PTR)  :  INTEGER; 

(*  Calculates  the  cost  of  the  fan  of  the  subpolygon  *) 

(*  bounded  above  by  HC  and  below  by  HR.  *) 

VAR  TEMPI,  TEMP 2  :  INTEGER; 


BEGIN 

TEMPI  ;=  HR" . S IDE_PRODUCT  -  HC* . S IDE_PRODUCT 

+  HC* . BASE_PRODUCT; 

WITH  HR"  DO 

IF  HEAD_SMALL 
THEN 
BEGIN 

IF  HEAD  =  HC" .HEAD 

THEN  TEMP 2  :=  HC" . BASE_PRODUCT 

ELSE  TEMP2  :=  CP[HEAD  +  lj  -  CP  [HEAD]  ; 

FAN_COST  :=  (TEMPI  -  TEMP 2 )  *  W [ HEAD] ; 

END 

ELSE 

BEGIN 

IF  TAIL  =  HC* . TAI L 

THEN  TEMP 2  ;=  HC* . BAS E_PRODUCT 

ELSE  TEMP2  :=  CP  [TAIL]  -  CP[TAIL-1J; 

FAN_COST  :=  (TEMPI  -  TEMP 2 )  *  W[TA1I,1; 

END; 

END;  {fan  cost} 


FUNCTION  SUPPORT  I  NG_WE  IGNT  (HR,  HC  :  L.IST_PTR)  :  INTEGER; 

(********************************************************) 
(*  Find  the  supporting  weight  of  the  subpolygon  bounded  *) 
(*  above  by  HC  and  below  by  HR.  *) 

^★*******************************************************) 
VAR  Y  :  INTEGER; 


BEGIN 

{calculate  the  denominator} 

Y  :  =  (HR*  .  SI  DE__PRODUCT  -  HR*  .  RAGK_PRODUCT) 

-  (IIC*. SIDE  PRODUCT  -  HO*  .  BASK  PRODUCT); 


{calculate  the  SUPPORT!  NG_WE  IGIIT } 

SU  Pl’O  RT I  NG_WE  IGHT  ;=  (HR*. COST  +  Y  -  1)  DIV  Y; 

{ceil i ng  func  t ion } 

END ;  {  s uppo r  t  i  ng__we  ig  h  t } 
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PROCEDURE  REMOVE  (VAR  S  :  r.lST  PTE;  MIN  :  INTEGER)  ; 

(*******  ***************  tlHUU'HHHHiHttinHH - •  -  -  j 

(*  Removes  al  l  the  arcs  in  S  whose  SU 1 ’_W LIGHTS  an*  equal  to  *  ) 
(*or larger  than  MIN  from  the  list.  *) 

^********************************************************J 

VAR  NOT  DONE  :  BOOLEAN; 


BEGIN 

NOT  DONE  :=  TRUE; 


WHILE  NOT_DONE  DO 
IF  S  =  NIL 

THEN  NOT_DONE  :=  FALSE 
ELSE 

IF  S' .SUP_WEIGHT  <  MIN 
THEN  NOT_DONE  :=  FALSE 
ELSE  S  :=  S' .NEXT; 

END;  {remove} 


PROCEDURE  SUB_PARTITION  (VAR  S  :  LIST_PTR;  MIN  :  INTEGER); 

(*  Finds  the  optimum  partition  of  the  subpol ygonbounded  *) 
(*  below  by  the  potential  h-arc  at  the  head  of  S.  The  *) 
(*  h-arcs  in  the  optimum  partition  of  the  subpolygon  *) 
(*  is  kept  in  a  list  with  S  pointing  to  the  head  of  *) 
( *  the  list.  * ) 

(********************************************************) 
VAR  TEMP  ;  INTEGER; 

TEM P_PTR  :  LIST_PTR; 

NOT_DONE  :  BOOLEAN; 

BEGIN 

IF  S'. NEXT  <>  NIL 
THEN 
BEGIN 

IF  S'  . HEAD_SMAI.L 
THEN  TEMP  :=  W [S'. HEAD] 

ELSE  TEMP  :=  W[S'.TAIL]; 

SUB_PARTITION  (S' .NEXT, TEMP) ;  {S'. NEXT  may  become 

Nil,  when  return 
from  S U B_P A RT I TON } 

END; 


IF  S'. NEXT  =  Nil, 

THEN  TEMP_PTR  ;=  LEAF  (S  is  the  last  arc  in  the  list] 

{ LEAF  is  a  psuedo  arc  with 
both  LEAF'. RASE  PRODUCT  and 
LEAF'  .SID E_ P RO I )UCT  equal  to  NIL} 
ELSE  TEMP  PTR  :=  S'.NEXT; 
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. COST  :=  FAN_COST(S ,TEMP_PTR) ; 

NOT  I  JONH  TRUE; 

WHILE  NOT_UONH  IX) 

BHGIN 

S' . SUP_WE IGHT  :=  SUPPORT I NG_WE IGHT ( S , THM P_PTR) ; 

IF  S“ .SUP_WEIGHT  >  =  MIN  (to  see  if  the  partition  is 

optimum  in  the  subpolygon} 

THEN 

BEGIN 

REMOVE  (S,MIN);  {delete  all  h-arcs  not  in  the 

optimum  partition  of  the 
subpolygon } 

NOT_DONE  :=  FALSE; 

END 

ELSE 

BEGIN 

IF  S~ .NEXT  <>  NIL 
THEN 

IF  S~  .NEXT~  .SUP_WETGIIT  <=  S~  . SUP_WE IGHT 

THEN  NO T  DONE  FALSE 

ELSE 

BEGIN  {condense  S' .NEXT  into  £} 

TEMP_PTR  :=  S~ .NEXT; 

S"'  .NEXT  :=  TEMP_PTR“  .NEXT; 

S''.COST  :=  S*  .COST  +  TEMP_PTR*  .COST; 
TEMP_PTR~ .NEXT  :=  S~ .DESCENDANT; 

S' . DESCENDANT  :=  TEMP_PTR; 

IF  S' .NEXT  =  NIL 
THEN  TEMP_PTR  LEAF 
ELSE  TEMP_PTR  :=  S~  .NEXT; 

END 

ELSE  NOT_DONE  :=  FALSE; 

END; 

END; 

END;  {sub_par tition} 


BEGIN 

S U B_ P A RT I T 1 0 N  (L  ,  W{  1  ]  )  ; 
END;  { mono_par tition} 
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PROCEDURE  WRITE  LIST  (L  :  LIST  PTR ;  MIN,  INDENT  :  INTEGER); 
(*************  *1F  *************  *'f  ****************************  ) 

( *  Displays  the  h-arcs  in  the  list  pointed  by  L.  *) 

|**********************************************ik*********ak*j 

VAR  TEMP  :  POS_I NTEGER ; 

BEGIN 

WHILE  L  <>  NIL  DO 
BEGIN 

IF  L*  .HEAD_SMALL 
THEN  TEMP  :=  L~.HEAD 
ELSE  TEMP  :=  L'.TAIL; 

IF  TEMP  <>  MIN 

THEN  WRITELN  ('  ' ; INDENT , lT . HEAD , '  ' : 3 , L* .TAIL) ; 

WRITE__LI ST  (L'  .  DESCENDANT  ,  TEMP  ,  I NDENT+  3  )  ; 

L  :=  L~ .NEXT; 

END; 

END;  {write  list} 


BEGIN  {main  program  begins  here} 
INITIALIZING; 

ON  E_SWEEP  (LIST); 
MONO_PARTITION  (LIST); 

IF  LIST  <>  Nil. 

THEN  WRITE_LIST  (LIST,  1,3) 

ELSE  WRITELN  (’  ’ :3,'NIL'); 

END.  {main  program} 


Appendix  II 


P  HOG  R A M  0  PT 1 M  U  M_P  A  RT I T I O  N_0  E_ A_G  E  N  K  R A  L_C  O  N  V  E  X_PO  L Y  GO  N  ; 

CONST  MAX_SIZE  =  127;  {the  maximum  number  of  vertices  in 

a  polygon  is  126} 

MAX  INT  =  32767;  {the  largest  integer  in  the  machine} 


TYPE  POS_I NTPGER 
LIST_PTR 
LIST  ELEMENT 


TREE_PTR 
TREE  ELEMENT 


0  ..  MAX_I NT ; 
“L1ST_EI.EMKNT 
PACKED  RECORD 
HEAD 
STAY 
TAIL 

HEAD_SMALL 

NEXT 

END; 


POS_INTEGER; 
BOOLEAN ; 
POS_INTEGKR; 
BOOLEAN; 

LIST  PTR 


~TREE_ELEMENT; 

PACKED  RECORD 

H EAD ,  TA I L  : 

HEAD__SMAI,I.  : 

SUP_WEIGHT, 

TREE_COST, 

TREE_BASE_PRODUCT, 

TRE E_S  1  f)E_PRODUCT  , 

LOCAL_COST, 

LOC  A I  j_  B  A  S  E_  PRODUCT  , 
LOCAL_S  I  DE_PROD(JCT : 
DESCENDANT,  NEXT, 
H_ARC,  V_ARC  : 

LIST_LINK  : 

DEPTH  : 


POS_I NTEGER ; 
BOOLEAN; 


INTEGER; 

TREE_PTR; 
LI ST_PTK ; 
INTEGER 


END; 


{V_ARC  and  H_ARC  are  used  in  two 
different  ways  :  (1)  they  are 

used  to  link  the  unprocessed  arcs 
together  to  form  an  arc-tree;  and 
(2)  they  are  used  as  the  left 
link  and  the  right  link  of  the 
processed  arcs  in  the  leftist 
tree  for  the  priority  queue.  } 


VAR  W,  CP 

1, 1ST  l  ,  LI  ST  2 
VJTREE,  H_TREE 
N 


ARRAY  ( 1 . .MAX_SIZK]  OE  INTEGER; 
r,TST_PTR; 

TRKK_PTR ; 

POS  INTEGER; 
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SEGMENT  PROCEDURE  INITIALIZING; 

|*****************************i**-A*********'******A********ij 

(*  Handles  the  inputs  and  ini  t  i  a  1  i  z  i  ng  all  the  cjlob.il  *) 

(*  variables.  *) 

(**********************************************************) 
VAR  I  :  INTEGER; 


BEGIN 

WRITET.N  ('a  linear  algorithm  to  find  all  the  h-arcs  in'); 
WRITELN  ('  the  optimum  partition  of  a  convex  polygon', 

'  (7/16/80)'); 

WRITELN; 

{obtain  the  inputs) 

WRITE  ('Please  enter  the  size  of  the  polygon  (between  3', 

’  and  '  ,MAX  SIZE-1  ,  ' )  ;  ' )  ; 

READI.N  (N); 

WRITELN; 

WRITELN  ('Now,  starting  from  the  smallest', 

'  vertex  and  in  the  '); 

WRITELN  ('  clockwise  direction,  ', 

'enter  the  weights  of  the  vertices:'); 
FOR  I  :=  1  TO  N  DO  READ  ( W [ I )  )  ; 

READLN ; 

WRITELN; 

(calculate  the  cumulative  adjacent 

products  around  the  polygon} 

CP [ 1 ]  :=  0; 

FOR  I  :=  2  TO  N  DO  CP  [II  :  =  CP  ( I  - 1  ]  +  W  [  I  - 1  ]  *  W  [  I J  ; 

{set  up  the  output  headings} 

WRITELN  ('the  potential  h-arcs  in  the  partitions  are  :  *); 
END;  {initializing} 
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mmmmm 


segment  procedure  one  sweep  (vak  l  :  list  ptr)  ; 

(  A  A  A  A  A  A  A  A  AAA*  A*  A  A  A  A  A  A'f  **************  A  A  A  A  A  T  A  **********  AAA***) 

(*  Sweep  the  polygon  once,  collects  all  potential  h-arcs,  *) 
(*  puts  them  in  a  list.  The  address  of  the  head  of  the  *) 
(*  list  is  stored  in  I..  *) 

^***********************************A**************A******* 

VAR  STACK  :  ARRAY  [ 1 . . MAX_S I ZE 1  OF 

POS_INTEGER; 

TOP_EI.EMENT ,  S ECON  D_E L EM F.NT  , 

CURRENT,  TOS  :  POS_I NTEGER ; 

P,  ARC  LIST  :  LIST  PTR ; 


PROCEDURE  PUSH  (C  :  INTEGER)  ; 


(*  Pushes  the  index  C  onto  the  stack  and  updates  the  *) 

(*  variables  TOS,  TOP  ELEMENT,  and  SECOND  ELEMENT.  *) 


BEGIN 

STACK [TOS]  :=  C; 

SECOND_ELF,MENT  :=  TO  P_E  L  EM  ENT ; 
TOP_ELEMENT  ;=  C; 

TOS  :=  TOS  -  1; 


END;  [push] 


PROCEDURE  POP_STACK 
( *  *  * 

(*  Pops  the  top  element  off  the  stack  and  updates  the 
(*  variables  TOS,  TOP_E  LEM  ENT ,  and  SECOND  ELEMENT’ . 

( 


*  *  *  *  *  J 
*) 
*) 
) 


BEG  1  N 

TOS  :=  TOS  +  1; 

TOP_ELEMENT  SECON  D_E  LEMENT ; 

S  ECOND_F.  LEM  F.NT  ;=  STACK  [TOS  +  2]; 

END;  lpop_Stack) 


BEGIN  [one  sweep  begins  here} 

{initialTze  the  local  variables) 
TO  P__  ELEMENT  ;=  0; 

SECOND  ELEMENT  ;  =  0; 

STACK [Nil]  :=  0; 

TOS  :=  N; 

ARC_L1ST  ;=  NIL; 

PUSH  (1); 

PUSH  (2) ; 

CURRENT  ;=  3; 
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{scan  through  the  polygon  in  the  clockwise  direction} 

WHILE  CURRENT  <  N  DO 

IK  (W [SECOND  ELEMENT)  <;  W[TOP  ELEMENT)  )  AND 

( W [ TO P  ELEMENT)  >  W  [CURRENT)  ) 


THEN 
BEGIN 
NEW ( P ) ; 

WITH  P“  DO 
BEGIN 

HEAD  :=  SECOND  ELEMENT; 

TAIL  :=  CURRENT; 

STAY  :=  FALSE; 

HEAD_SMALL  :  =  W {HEAD]  <=  WlTAIL}; 
NEXT  :=  ARC_LIST; 

END; 

ARC  LIST  :=  P; 


POP_STACK; 

IF  lOS  >=  (N  —  1 )  {there  are  less  than 

2  elements  on  the  stack} 

THEN 

BEGIN 

PUSH  (CURRENT) ; 

CURRENT  :=  CURRENT  +  1; 

END; 

END 

ELSE 

BEGIN 

PUSH  (CURRENT) ; 

CURRENT  :=  CURRENT  +  1; 

END; 

WHILE  (TOS  <=  (N- 3 )  ) 

AND  ( W  [  S  ECO N D_F.  LEM E NT]  <=  W  [  T O P_K I . K M E N T 1  ) 

AND  ( W [ TO P_E LEMEN T)  >  W[N])  DO 

BEGIN 

NEW(P)  ; 

WITH  P'  DO 
BEGIN 

HEAD  :=  S  ECOND_E  LEM ENT ; 

TAIL  :=  N; 

STAY  :=  FALSE; 

HEAD_SMALI.  :=  W(HEAD)  WlTAIL); 

NEXT  :=  ARC_I.IST; 

END; 

ARC_LIST  :=  P.; 

POP  STACK; 

END; 


L  ;=  A  RC_L I  ST ; 
END;  {one_sweep} 
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segment  procedure  RUILD_TREE  (VAR  I.  :  LIST  PTK; 

VAR  VT,  HT  :  TREE_PTR;  FIRST,  LAST,  mTn  :  POS  INTEGER); 
(***** ******************************************* *V***  A****) 

(*  Traces  all  the  arcs  in  the  list  pointed  by  L  and  *) 

(*  build  an  arc-tree  with  the  root  pointed  by  T.  *) 

(**********************************************************) 
VAR  NOTJDONE  ;  BOOLEAN; 

P  :  TREE_PTk; 

Q  :  LI ST_PTR ; 

BEGIN 

NOT_DONE  :=  TRUE; 

VT  :=  NIL; 

HT  ;=  NIL; 

WHILE  NOT_DONE  DO 
IF  L  =  NIL 

THEN  NOT_DONE  :=  FALSE 
ELSE 

IF  (I, ".HEAD  <  FIRST)  OR  (L~.TAIL  >  LAST) 

THEN  NOT_DONE  :=  FALSE 
ELSE 
BEGIN 

Q  ;=  L" .NEXT; 

IF  L* . HEAD  <>  1 
THEN 
BEGIN 

NEW  (P); 

WITH  P~  DO 
BEGIN 

HEAD  :=  L~ .HEAD; 

TAIL  :=  L* . TAIL; 

HEAD_SNALL  :=  L~ . HEAD_SMALL; 

DESCENDANT  :=  NIL; 

DEPTH  :=  1; 

LIST_LINK  ;=  L; 

{ LOCAL_COST  ,  LOCAI,_BASE_PRODUCT , 

LOCAL  SID E_P RODUCT ,  TREE_COST, 

T RE E_B  AS E_P ROD U C T  ,  TREE_S  I  l)E_P RODUCT  , 
H__ARC,  and  V  ARC  are  undefined  at  this 
point*"  } 

IF  (HEAD_SMAI.L  AND  {HEAD  =  MIN)  )  OR 

(NOT  HEAD_SMALL  AND  (TAIL  =  MIN)) 

THEN 

BEGIN 

NEXT  VT; 

VT  ;=  P; 

END 

ELSE 

BEGIN 

NEXT  :=■  II T; 

HT  :=  P; 

END; 
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IF  HEAD  SMALL 
THEN  lUlTl,!)  THEE  (0,V  ARC,  II  ARC, 

"HEAD, TAIL, HEAD) 
ELSE  BUI  l.D_TREE  (Q ,  V_ARC  ,  H_ARC  , 

HEAD, TAIL, TAIL) ; 

{note  that  there  will  be  at  most  one  ar^. 
in  the  V_ARC  list  but  may  be  several  arcs 
in  the  H__ARC  list  } 

END; 

END; 

L  :=  Q; 

END; 

END;  {build  tree} 


SEGMENT  PROCEDURE  POLY  PARTITION  (VAR  T  :  TREE  PTR) ; 


**********; 


(*  To  find  all  the  h-arcs  that  are  present  in  the  optimum  *) 
(*  partition  of  the  polygon  and  returns  them  in  the  arc-  *) 
(*  tree  pointed  by  T.  *) 


PROCEDURE  FAN_COST  (T  :  T'RE E_PT R 

(*  To  find  the  cost  of  the  fan  of  the  subpolygon  bounded*) 
(*  below  by  the  arc  pointed  by  T  and  above  by  the  arcs  *) 
(*  pointed  by  T~.H__ARCs  and  T~.V_arcs.  *) 

^★********************A***************A******************^ 

VAR  X  ;  POS_I NTEGE R ; 

Y,  SI,  S2  :  INTEGER; 


BEGIN 

WITH  T“  DO 
BEGIN 

IF  HEAD  SMALL 
THEN 
BEGIN 

IF  V_ARC  =  NIL 
THEN 
BEGIN 

X  :=  HEAD  +  1; 

SI  :=  CP{X)  -  CP [ HEAD) ; 

END 

ELSE 

BEGIN 

X  V_ARC".TAIL; 

SI  :=  V_ARC* . TREE_BASE_PRODUCT ; 
END; 
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S 2  :=  (CP[TAJI,|  -  CP[X)  )  ; 

Y  :  -  W[Ili:ADl  ; 

END 

ELSE 

BEGIN 

IP  V_ARC  =  NIL 
THEN 
BEGIN 

X  :=  TAIL  -  1; 

SI  :=  CP  [  TAI  Ll  -  CP[X]  ; 

END 

ELSE 

BEGIN 

X  :=  V_ARC" .HEAD; 

SI  :=  V_ARC~ . TREE_BASE_PRODUCT ; 

END; 

S  2  :=  (C P  ( X ]  -  CP[HEAD]  )  ; 

Y  :=  W [ TA I L ] ; 

END; 

IP  H_A RC  <>  NIL 

THEN  S 2  :=  S2  -  H_ARC~  .TREE_S1DE  PRODUCT 

+  H_ARC’r.TREE_BASE_PRODUCT; 
[all  the  SI DE_PRODUCTs  and  the  BASEJPRODUCTs  are 
added  together  and  stored  in  the  root  of  the 
leftist  tree  pointed  by  H_ARC  } 

LOCAL_COST  :=  S2  *  Y; 

LOCAL_SIDE_PRODUCT  :=  Si  +  S2; 

LOCAL_BASE_PRODUCT  :=  W[HEAD)  *  W[TAIL]; 

END; 

END;  [fan  cost} 


PROCEDURE  SUPPORTING_WK IGHT  (T  :  TR E E_  P T R )  ; 

^  ********************************************************  j 

(*  To  find  the  supporting  weight  of  the  arc  pointed  by  T* ) 
(*  with  respect  to  the  suhpolygon  bounded  below  by  the  *) 
(*  arc  pointed  by  T  and  above  by  the  arcs  pointed  by  *) 
(*  the  T^.H^ARC  and  T^.V^ARC.  *) 

VAR  D  :  INTEGER; 

BEG  I  N 

WITH  T~  DO 
BEGIN 

D  ;=  (LOCA I. _ S  1  DE_PROUUCT  -  LOC AL_BASE_PUODUCT)  ; 

SUP_WEIGHT  :=  ( LOCAL_COST  +  D  -  1)  D1V  D; 

(ceiling  function} 

END; 

END;  (suppor  ting_we  ight } 
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FUNCTION  MERGE  (Tl,  T2  :  TUFF  PTH)  :  TREEPTK; 

^***** *************  *********  *T***********f***************j 

(*  Merqes  two  leftist  trees  inLo  one  and  returns  it  in  *) 
(*  MERCK .  *) 

^********************************************************) 
VAR  TKMP_PT1<  :  TREK  P'l'R; 

TEMP_COST,  TKMP_BASK_PRODUCT, 

TEMP  SIDE  PRODUCT  :  INTEGER; 


BEGIN 

IF  T2  =  NIL 
THEN  MERGE  :=  Tl 
ELSE 

IF  Tl  =  NIL 
THEN  MERGE  :=  T2 
ELSE 
BEGIN 

TEMP_COST  :=  Tl" . T  RE  E_COST  +  T2* . TREE_COST ; 
TEMP_SIDE_PRODUCT  :=  T 1  *  . TRE  E_S I D  E_P  RO  DUCT 

+  T2* . TREE_S IDE_PRODUCT; 
TEMP_BASE_PRODUCT  :=  Tl* .TREE  BASR_PRODUCT 

+  Tl * . T  R  E E_B  A  S  E_P  ROD  U  C  T ; 
IF  Tl* .SUP_WEIGHT  <  T2* . SUP_WE 1GHT 
THEN 
BEGIN 

T EM P_PT R  Tl  ; 

Tl  :=  T2; 

T2  :=  T  EM  P_PT R ; 

END ; 

WITH  Tl*  DO 
BEGIN 

H_ARC  :=  MERGE  (H__ARC,T2); 

{ H_ARC  never  equals  NIL  at  this  point} 

IF  V_ARC  =  NIL 
THEN 
BEGIN 

V  ARC  H  ARC; 

H~ARC  NTi. 

END 

ELSE 

BEGIN 

IF  V_A RC* . DEPTH  <  H  ARC*  . DEPTH 
THEN 
BEGIN 

TEMP__PTR  :=-  V  A  RC  ; 

V  ARC  :  =  11  ARC; 
h”ARC  TEMP  P'l'R; 

END; 

DEPTH  ;=  H_ARC*. DEPTH  +  1; 

END; 

TREE  COST  TEMP  COST; 
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# 


TREE_S  I  DE_PRODUCT  :=  TEMP_S 1 DE_PROPUCT; 
TREE  HAS E_PRODUCT  :=  TEMP  HAM-:  PRODUCT ; 
END; 

MERGE  :=  Tl; 

END; 

END;  {merge} 


FUNCTION  CONDENSE  (T  :  TRE E_PTR ;  MIN  :  INTEGER)  :  BOOLEAN; 

(*  CONDENSE  is  set  to  false  if  T  =  NIL  or  T~ . SUP_WE IGHT  *) 
(*  <=  MIN  *) 

|********************************************************-j 

BEGIN 

IF  T  =  NIL 

THEN  CONDENSE  :=  FALSE 

ELSE  CONDENSE  ;=  T“ .SUP_WE IGHT  >  MIN; 

END;  {condense} 


PROCEDURE  COMBINE  (VAR  T  :  TREE_PTR;  V_FI.AG  :  BOOLEAN)  ; 

^********************************************************) 
(*  If  V.FLAG,  it  combines  the  arc  pointed  by  T~.V_ARC  *) 
(*  into  the  arc  pointed  by  T,  else  it  combines  the  arc  *) 

(*  pointed  by  T".H_ARC  into  the  arc  pointed  by  T.  In  *) 

(*  either  case,  the  arc  to  be  combined  is  deleted  from  *) 

(*  the  corresponding  leftist  tree  and  put  into  the  *) 

(*  DESCENDANT  list  of  the  parent.  *) 

^■k^-kic-k-k-k'k-k-k-k^ic-k-k-kic^c-k'k'k-k-k-k-kic-k-k'kie'k'k'kic-k-kif'k'k-k'kit-kicicicic'kic'k^icAit-kif'j 

VAR  TEMP  PTR  :  TREE  PTR ; 


BEGIN 

IF  V_FLAG 
THEN 
BEGIN 


T EM P_PT R  :  =  T*  .  V_A RC ; 

T~.V_ARC  :=  MERGE  (TEMP  PTR"  .  V_ARC  ,  TKMP_PTR*  .  H  ARC); 
END 
ELSE 
BEGIN 


TFM P_PTR  :=  T".H_ARC; 

T "  .11  ARC  :=  MERGE  (TEMP  PTR‘.V  ARC, TEMP  PTR*. II  ARC); 
END;  “  " 


TEMP_PTR"  .  V_ARC  NIL; 

TEMP_PTR" ,H_ARC  :=  NIL; 

TEMP_PTR* .NEXT  T" .DESCENDANT; 

T" . DESCENDANT  :=  TEMP_PTK; 

T"  .I.OCAL_COST  T". LOCAL  COST  I  TEMP  PTK~  .  LOCAL  COST; 

T"  .  LOCAL_SIDE_PRODUCT  i';"  .  LOCAI._S I  DE_PRODUCT 
t  TEMP_PTR"  . LOCATES  1 DE_PRODUCT 

-  TEMP_PTR"  .  LOCAL _RASE_PROI>UCT ; 

END;  {combine} 
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PROCKIHIKR  RKMOVK  (VAR  T  :  TKRR_PTU;  MIN  :  1NTKGKR)  ; 

|***************H**H********H******HH****HxH**H**j 

(*  Removes  all  the  arcs  in  the  leftist  tree  pointed  by  T*  ) 
(*  who.se  SUP_WE  I  GUTS  are  larger  than  or  equal  to  MIN.  *) 

^*Ar**A**AAAAA*A-*AA*******AAAAAA*AAAAAA*x**A*AAAAAA*AAAr*A,A^ 

VAR  NOT_DONE  :  BOOLEAN ; 

BEGIN 

NOT_DONK  :---  TRUE; 

WHILE  NOT_DONE  DO 
IF  T  =  NIL 

THEN  NOT_DONE  :=  FALSE 
ELSE 

IF  T*  .SUP  WEIGHT  <  MIN 
THEN  NOT_DONE  :  =  FALSE 
ELSE  T  :=  MERGE  (T~ .  V_ARC,T~ . H_ARC) ; 

END;  {remove} 


PROCEDURE  SUB  PART  1'ION  {VAR  TREE 


THE E_PTR ; 

MIN  :  INTEGER)  ; 

^**A*A*****AA**Ar****Ar*A:AAAAAAAAAAAA*AAAAAAAAA*A**AAA*A:A:A:A;J 

(*  To  find  the  optimum  partition  of  the  subpolygon 
(*  bounded  below  by  the  root  of  the  arc- tree  pointed 
(*  by  T. 

|AA*AAAAAAA*AAAAAAAAAA*A*A*A***A  aa  *********** 

VAR  T,  R,  P,  TEM P_PTR  :  TREE_PTR; 

TEMP  ;  INTEGER; 

NOT  DONE,  FLAG  :  BOOLEAN; 


*) 

*) 

*) 

) 


BEGIN 

T  :=  TREE; 

R  :=  NIL; 

WHILE  TO  NIL  DO 
BEGIN 

P  :  =  T~  .NEXT; 

T* . NEXT  NIL; 

IF  T~ .HEAD_SMALL 
THEN  TEMP  :=  W  [T*  .  HEAD} 

ELSE  TEMP  :=  W[T* .TAIL) ; 

IF  T~ . H_ARC  <>  NIL 

THEN  SUBPARTITION  (T*  .  H_ARC  ,  TEMP)  ; 

{when  return,  all  the  h-ares  in  the  subpolygon 
will  be  put  in  a  priority  queue  } 

IF  T‘.V_ARC  <>  NIL 

THEN  S U B_P A RT I T ] O N  (T~ . V_AKC , TEMP)  ; 

{there  should  be  at  most  .1  v-arc,  i.e. 

T~ .V_ARC“ .NEXT  -  NIL,  when  return,  all  the 
h-ares  in  the  subpolygon  will  be  put  in  a 
priority  queue  } 
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{calculate  the  cost  of  the  fan  of 
bounded  below  by  the  arc  po  ntr-d 
the  v-arcs  and  h-urcs  of  t 
FAN_COST (T) ; 

NOT  DONE  :=  TRUE; 

FLAG  :=  TRUE; 

WHILE  NOT_DONE  DO 
BEGIN 

{calculate  the  supporting  weight 

of  the  arc  pointed  by  T} 

SU PPORTING_WE IGHT  (T)  ; 

IF  T" ,SUP_WE1GHT  >=M1N  {to  see  if  the  partition 

is  optimum  in  the 
subpolygon  } 

THEN 

BEGIN 

REMOVE  (T,MIN);  {delete  all  h-arcs  not 

in  the  optimum  partition 
of  the  subpolygon  } 

NOT_DONE  :=  FALSE; 

FLAG  :=  FALSE; 

END 

ELSE 

IF  CONDENSE  (T~  .  V__ARC  , T~  .  SUP_WE IGHT) 

THEN  COMBINE  (T ,  TRUE) 

ELSE 

IF  CONDENSE  |f . H_ARC , T* . SUP_KEIGHT) 

THEN  COMBINE  (T,  FALSE) 

ELSE  NOT  DONE  FALSE; 

END; 

{maintain  the  leftist  tree  structure] 

IF  FLAG 
THEN 
BEGIN 

T*  .TREE_COST  :=  T* . LOCALJTOST; 

T- .TRBE_S ID BYPRODUCT  . LOCAL_SIDE_PRODUCT; 

T~  .TREE_BASE_PRODUCT  T*  .  LOCAL  BASE_P RODUCT  ; 

IF  T~ . V_ARC  <>  NIL 
THEN 
BEGIN 

Ta ,TREE_COST  :=  T* .TREE_COST 

+  Ta.V  ARC" .TREE_COST; 
T"  .TREE_SIl)E  PRODUCT  TA.TRF.L  SI  I  )!■:_  PRODUCT 
+  Ta‘."V_ARC' .TREE  SIDE  PRODUCT 

-  T"  .V_ARC’'.TRKE_BASF._PRODUCT; 

END; 


the  subpojyyon 
by  T  and  above  by 

} 
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r  -  - 

: 

IP  T“.H_ARC  <>  Nil. 

THEN 

k  BMC  IN 

T~  .TRl-:i:_COST  :  =  T*.TREK  COST 

+  "'t~.!1  ARC"*  .TRKE_COST; 
t“  .TKEi:_sii)i:_i>uoi)(jCT  t"  .  tki;f:_si  I)i;_phoouct 
+  T*.H_ARC~ .TREE  S  I  DK_PRODUCT 

-  T'  .H_ARC'.TREE_BASE_PRODUCT; 

END; 

A- 

IF  T  . V_A RC  <>  NIL 
THEN 

IF  T~  .  H_A  RC  <>  NIL 
THEN 
BEGIN 

E  IF  T* . V_ARC” . DEPTH  <  T~ . H_ARC~ . DEPTH 

THEN 

BEGIN 

TEMP_PTR  :=  T*.V_ARC; 

T".V_ARC  :=  T“.H_ARC; 

T"  .H_ARC  :  =  TEMP_PTR; 

END; 

T~ . DEPTH  ;=  T~ . H_ARC~ . DEPTH  +  1; 

,  END 

ELSE 

ELSE 

IF  T~.H_ARC  <>  NIL 
THEN 
BEGIN 

T“.V_ARC  :=  T* . H_ARC ; 

T".H_ARC  :=  NIL; 

END; 

END; 


R  :=  MERGE  { R  ,  T)  ; 
T  :=  P; 

END; 

TREE  :=  R; 

END;  (subpartition) 


BEGIN  (poly  partition  begins  here) 
SUB  PARTITION  (T,W(  1  ]  )  ; 

END;  Tt/Uy_parti  tionj 


I 

I 

I 
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PROCEDURE  MARK  1. 1ST  (T  :  TREE  PTR)  ; 

(************  *T  *************  *T  **************************  ***j 

(*  Traverses  the  tree  pointed  by  T  preordci  ly,  finds  out  *) 

(*  all  the  potential  h-arcs  which  are  present  in  tlie  *) 

(*  optimum  partition  of  the  polygon  and  marks  the  *) 


BEGIN 


WHILE  T  <>  NIL  DO 
BEGIN 

T“ . LIST_LINK" .STAY  :=  TRUE; 
MARK_LIST  (T" . DESCENDANT) ; 
MARK_LIST  (T" .V_ARC) ; 
MARK_LIST  (T".H_ARC); 

T  :=  T~ .NEXT; 

END; 

END;  {mark  list} 


PROCEDURE  WRITE_LIST  (VAR  L  :  LISTJPTR; 

FIRST,  LAST,  MIN,  INDENT  :  INTEGER); 

(**********************************************************) 
(*  Displays  the  h-arcs  in  the  list  pointed  by  L.  *) 

(**********************************************************} 
VAR  TEMP  :  POS_INTEGER; 

NOT_DONE  :  BOOLEAN; 

BEGIN 

NOT_DONE  :=  TRUE; 

WHILE  NOT_DONE  DO 
IF  L  =  NIL 

THEN  NOT_DONF.  :=  FALSE 
ELSE 

IF  ( L~  .HEAD  <  FIRST)  OR  (L".TAIL  >  LAST) 

THEN  NOT_DONE  :=  FALSE 
ELSE 
BEGTN 

IF  L"  .STAY 
THEN 
BEGIN 

IF  L"  . HEAD_SMALL 
THEN  TEMP  :=  I."  .HEAD 
ELSE  TEMP  ;=  L~ .TAIL; 

IF  TEMP  <>  MIN 
THEN 
BEGIN 

WRITE LN  (’  ’ ; INDENT, 


L"  .  HEAD  , '  f:  3,1.".  TAIL) 


WRITE_LIST  (L  .NEXT,:.  .HEAD, 

L~ .TAIL, TEMP, INDENT*-!) 

END; 

END; 

L  :=  L~ .NEXT; 

END; 

END;  (writ e_ list} 


(*  main  program  begins  here. 

(**************************************** 

BEGIN 

INITIALIZING; 

ON  E_SWEEP  ( LISTl )  ; 

LIST2  :=  LISTl; 

BUILD_TRF,E  ( L I  ST  2  ,  V_T  REE  ,  H_T  REE  ,  1 ,  N  ,  1 )  ; 
POLYJPARTITION  (H_TREE)  ; 

IF  H_TREK  =  NIL 

THEN  WRITF.LN  ('  ':3,'NIL') 

ELSE 

BEGIN 

MARK_1. 1ST  (H_TREE)  ; 

WRITE_LIST  (LISTl ,1,N,1 ,3) ; 

END; 

END.  (main  program} 


{V  TREE  =  NIL} 
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